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ABSTRACT 

We present single-epoch, post-maximum spectropolarimetry of four Type 
la supernovae (SNe la) that span a range of spectral and photometric proper- 
ties: SN 2002bf and SN 2004dt exhibit unusually high-velocity (HV) absorp- 
tion lines. SN 1997dt is probably somewhat subluminous, and SN 2003du is 
slightly overluminous. We detect polarization modulations across strong lines 
in all four objects, demonstrating that all are intrinsically polarized. How- 
ever, the nature and degree of the polarization varies considerably. Including 
all SNe la studied thus far, the following order emerges in terms of increasing 
strength of line-polarization features: ordinary/over luminous < subluminous < 
HV SNe la, with the strength of the line-polarization features increasing from 
0.2% in the slightly overluminous SN 2003du to 2% in both of the HV SNe la 
in our study. The most convincing explanation for the line polarization of all 
objects is partial obscuration of the photosphere by clumps of intermediate- 
mass elements (IMEs) forged in the explosion; the polarization characteristics 
of the HV SNe la in particular effectively rule out a simple ellipsoidal aspheric- 
ity as the root cause of their line polarization. That SNe la are separable into 
different groups based on their spectropolarimetric characteristics may help nar- 
row down progenitor possibilities and/or explosion physics. In particular, the 
recently proposed gravitationally confined detonation model may provide an 
attractive explanation for many of the observed polarization characteristics of 
HV SNe la. 

Subject headings: supernovae: individual (SN 1997dt, SN 2002bf, SN 2003du, 
SN 2004dt) — techniques: polarimetric 
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1. Introduction 

The substantial homogeneity of Type la supernovae (SNe la), together with the fact 
that the peak luminosity of individual objects can be observationally determined through 
comparisons with well-calibrated samples of nearby SNe la, has propelled them to "gold 
standard" status among extragalactic distance indicators. Precise distance measurements 
out to 2; 1.7 have been made (Riess et al. 2001, 2004), revealing the surprising cosmo- 
logical result that the expansion rate of the universe is currently accelerating (Riess et al. 
1998a; Perlmutter et al. 1999); see Filippenko (2004, 2005) for extensive reviews. Although 
the progenitor systems have not yet been conclusively identified, the general consensus is 
that SNe la arise from carbon-oxygen white dwarfs (CO WDs) that accrete matter through 
some mechanism until they achieve a density of ~ 3 x 10^ g cm~^ in their centers, leading 
to a runaway thermonuclear reaction that incinerates the star (Woosley & Weaver 1986). 
This occurs, coincidentally, when a CO WD's mass is nearly the Chandrasekhar hmit of 
~1.4 Mq. Simulations of exploding CO WDs are able to reproduce the main spectral 
and photometric characteristics of SNe la (e.g., Leibundgut 2000, and references therein), 
lending theoretical support to this scenario. 

However, many questions remain concerning the SN la progenitors and explosion mech- 
anism. What is the nature of the "donor" that is responsible for providing the material that 
the WD accretes? Is there more than one channel by which the accretion can take place? 
Where docs the thermonuclear nmaway begin inside the WD, and how docs it propagate 
throughout the star? As with so many things, the "devil is in the details," and there is a 
growing sentiment that the answers to these fundamental questions may come from careful 
study of the differences seen among SNe la, rather than from the similarities alone. 

Spectropolarimetry offers the only direct probe of early-time SN geometry, and thus 
is an important diagnostic tool for discriminating among SN la progenitor systems and 
theories of the explosion physics. The essential idea is this: A hot, young SN atmosphere is 
dominated by electron scattering, which by its nature is highly polarizing. For an unresolved 
source that has a spherical distribution of scattering electrons, the directional components 
of the electric vectors of the scattered photons cancel exactly, yielding zero net linear 
polarization. Any asymmetry in the distribution of the scattering electrons, or of absorbing 
material overlying the electron- scattering atmosphere, results in incomplete cancellation, 
and produces a net polarization (see, e.g., Leonard & Filippenko 2005). 

Initial broad-band polarimetry studies found SNe la to possess zero or, at most, very 
weak (< 0.2%) intrinsic polarization (Wang et al. 1996), suggesting a high degree of sym- 
metry for their scattering atmospheres. More recent spectropolarimetric studies capable 
of resolving individual line features, however, are revealing a complex picture, with both 
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continuum and line polarization now convincingly established for at least a subset of the 
SN la population. 

To date, three SNe la have been examined in detail with spectropolarimetry at early 
times: SN 1999by (Howell et al. 2001), SN 2001el (Kasen et al. 2003; Wang et al. 2003), 
and, most recently, SN 2004dt (Wang et al. 2005a). In this paper we present single-epoch 
spectropolarimetry of four SNe la: SN 1997dt, SN 2002bf, SN 2003du, and SN 2004dt, 
obtained about 21, 3, 18, and 4 days (respectively) after maximum light. 

A particular motivation for this multi-object investigation is to attempt to link spectral 
and photometric peculiarities of individual SNe la with their spectropolarimetric character- 
istics. Accordingly, the four events span a range of properties: SN 1997dt is likely somewhat 
subluminous, SN 2002bf and SN 2004dt exhibit unusually high-velocity absorption lines (in 
the case of SN 2002bf, the highest ever seen in an SN la for the epochs considered), and 
SN 2003du is slightly overluminous. 

This paper is organized as follows. We briefly review the present photometric, spec- 
troscopic, and spectropolarimetric state of knowledge of SNe la in § 2, focusing particular 
attention on objects sharing characteristics with the SNe la included in our spectropo- 
larimetric survey. We describe and present the spectropolarimetry in § 3; to assist in the 
classification, we also include optical photometry and additional spectroscopy for two of 
the events, SN 2002bf and SN 2003du. We analyze the data in § 4, and present our conclu- 
sions in § 5. Note that preliminary discussions of the spectropolarimetry of SN 1997dt and 
SN 2002bf have been given by Leonard et al. (2000b) and Fihppenko & Leonard (2004), 
respectively. 

2. Background 

2.1. Photometric Properties of SNe la 

Type la SNe typically rise to peak 5-band brightness in about 20 days, decline by 
about 3 mag over the next 35 days, and then settle into a nearly constant descent of 
~1.55 mag (100 day)~^ for the next year. However, it is now confirmed beyond doubt 
that deviations exist from the central trend; see, for example, PhiUips et al. (1999). Some 
events rise and fall more slowly, producing "broader" light curves (e.g., SN 1991T), whereas 
others rise and fall more quickly, yielding "narrower" hght curves (e.g., SN 1991bg). The 
width of the light curves near peak correlates strongly with luminosity, in the sense that 
broader light curves generally indicate intrinsically brighter objects (but see Jha et al. 2005, 
for exceptions), with a total spread of about a factor of sixteen in absolute peak 5-band 
brightness among the population (Altavilla et al. 2004). 

The ability to individually determine the luminosity of SNe la through examination of 
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their light-curve shapes has been the primary driver for cosmological apphcations. There 
are several photometric calibration techniques in common use that correlate the observed 
light-curve shape, or supernova color, with luminosity (see Wang et al. 2005b, and references 
therein) . The simplest calibration technique involves measuring the decline in B during the 
first 15 days after maximum 5-band brightness, Amis (5). A typical value is Ami5(S) pa 
1.1 mag, which corresponds to Mb — —19.26 ± 0.05 mag according to the calibrations of 
Hamuy et al. (1996). Under luminous, SN 1991bg-like objects yield values up to Ami5{B) — 
1.94 mag, while overluminous, SN 1991T-hke objects go as low as Amir,(B) — 0.81 mag 
(Altavilla et al. 2004). There is a continuum of values in between the two extremes. 

For the two SNe with optical photometry presented in this paper, we shall use the "mul- 
ticolor light-curve shape" method (MLCS; e.g., Riess et al. 1996), which has been revised 
(and hereafter referred to as MLCS2k2) by Jha (2002) and S. Jha et al. (in preparation) 
to include f/-band light curves from Jha et al. (2005), a more self-consistent treatment of 
extinction, and an improved determination of the unreddened SN la color. This method is 
based on determining A (a dimensionless number related to magnitude that parameterizes 
the light-curve shape), to (the time of i?-band maximum light), Hq (the distance modu- 
lus for an assumed value of i?o, here taken to be 65 km s ~^ Mpc~^), and Ay (the visual 
extinction at to). 

The mechanism that produces the dispersion in SN la luminosity is not known; present 
speculations range from different progenitor systems to differences in the explosion mech- 
anism and flame propagation (e.g., deflagration, detonation, delayed-detonation, or, most 
recently, gravitationally conflned detonation; see, respectively, Nomoto et al. 1984; Arnett 
1969; Khokhlov 1991; Plewa et al. 2004). Some have further speculated that global asym- 
metries in the expanding ejecta may result in viewing-angle dependent luminosity (Wang 
et al. 2003). One observational fact that all theories must confront is that the most lumi- 
nous SNe la have thus far only been seen in late-type galaxies; the faintest objects tend to 
prefer elliptical galaxies, but have also been found in spirals (e.g., Howell 2001; Benetti 
et al. 2005, and references therein). 

2.2. Spectroscopic Properties of SNe la 

Early-time spectra of SNe la (see Filippenko 1997 for a review) typically exhibit lines of 
intermediate-mass elements (IMEs), such as magnesium, silicon, sulfur, and calcium, with 
some contribution from iron-peak elements, especially at near-ultraviolet wavelengths. As 
time progresses, and the photosphere recedes deeper into the ejecta, lines of Fe come to 
dominate the spectrum. This spectral evolution suggests a burning front that incinerates 
some of the progenitor's carbon and oxygen all the way to the iron peak deep inside the 
ejecta, but then leaves the outer layers only partially burned. It is interesting to note that 
overluminous SNe la show enhanced Fe features at early times, whereas subluminous ones 
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show weak early-time Fe features. 

At early times (e.g., —20 d < i < 20 d from the date of maximum B brightness, -Bmax), 
the spectra of typical SNe la evolve rapidly and with such uniformity that it is possible to 
determine the age of an event relative to the date of -Bmax to within ~2 days from a single 
spectrum alone (Riess et al. 1997; Foley et al. 2005). The spectroscopic peculiarities of 
subluminous and overhiminous events currently preclude their "spectral feature ages" from 
being derived accurately through comparison with average SN la spectra. For instance, 
pre-maximum spectra of overluminous events lack the strong Si II A6355 absorption that 
is so prominent in normal and underluminous SN la spectra (Filippenko et al. 1992b). 

Another point of distinction among SN la spectra comes from the blueshifts of the 
spectral lines. Increasing attention is being paid to a small but growing group of "high- 
velocity" (HV) SNe la, whose spectra around maximum light are characterized by unusually 
broad and highly blueshifted absorption troughs in many line features, indicating optically 
thick ejecta moving about 4000-5000 km s~^ faster than is typically seen for SNe la (Branch 
1987; Benetti et al. 2005). Well-studied examples of HV SNe la include SN 1983G (Branch 
& van den Bergh 1993, and references therein), SN 1984A (Branch 1987; Barbon et al. 
1989), SN 1997bp (Anupama 1997), SN 1997bq (Lentz et al. 2001, and references therein), 
SN 2002bo (Benetti et al. 2004), SN 2002dj (Benetti et al. 2004), and SN 2004dt (Wang 
et al. 2005a). Wang et al. (2005a) note that in a spectrum of one HV SN la, SN 2004dt, 
some line features do not possess abnormally high velocity, such as those identified with 
S II. Since two of our objects, SN 2002bf and SN 2004dt, are HV SNe la, we briefly review 
the sahent features that are known about this class of objects. 

The most thorough and recent study of the diversity of SN la expansion velocities is 
that of Benetti et al. (2005), who apply a statistical treatment to data from 26 SNe la and 
flnd that HV SNe la indeed make up a kinematically distinct group. Its members have 
normal peak luminosity (based on their Ami5{B) values: 1.09 < Ami5{B) < 1.37 mag), 
reside in all types of host galaxies (e.g., both ellipticals and spirals, but with a preference for 
later types), and have very strong Si II features. They are further distinguished by a large 
temporal velocity gradient, i; > 70 km s^^ day^^, where v is defined to be the average daily 
rate of decrease of the expansion velocity between maximum light and the time the Si II 
A6355 feature disappears. HV SNe la appear to have similar photometric characteristics 
to ordinary Type la events, although there is some indication of subtle differences in their 
color evolution (Benetti et al. 2004). 

Whether HV SNe la represent the extreme end of a continuum of more typical SNe la, 
or require a different explosion mechanism, progenitor system, or explosion physics, remains 
controversial. What is certain is that any model must have significant optical depth in the 
IMEs at high velocities at early times. Branch & van den Bergh (1993) originally proposed 
that HV SNe la may simply result from more energetic explosions. However, Benetti et al. 
(2004) show that the Fe nebular lines have velocities comparable to those of normal SNe la. 
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which tends to weaken the argument for higher overall kinetic energy. There is also evidence 
from early-time spectra of SN 2002bo that HV SNe la possess very little unburned carbon 
in their outer layers, in contrast with more typical SNe la like SN 1994D (Benetti et al. 
2004). The high velocities of the IMEs coupled with the lack of primordial carbon may 
indicate that burning to Si penetrates to much higher layers in HV SNe la than it does in 
more normal events. This is a feature of some delayed-detonation models, such as those 
studied by Lentz et al. (2001). However, Wang et al. (2005a) argue that a strong detonation 
wave is unlikely to generate the clumpy and asymmetrically distributed silicon layer that 
is inferred from pre-maximum spectropolarimetric observations of SN 2004dt. As another 
alternative, Benetti et al. (2004) propose that IMEs produced at deeper layers may simply 
be more efficiently mixed outward in HV SNe la than is typical. A prediction of this 
mechanism is that while the IME products of explosive carbon burning will be mixed in 
with the primordial C and O in the outer layers, the C and O should also be mixed inward 
within these events, and exist at lower velocities than are normally seen. 

It is appropriate at this point to mention the recently proposed gravitationally confined 
detonation (GCD) model of Plewa et al. (2004), which stands as an intriguing alternative 
to the standard deflagration or delayed-detonation scenarios. The mechanism involves the 
slightly off-center ignition of a deflagration that produces a buoyancy- driven bubble of 
material that reaches the stellar surface at supersonic speeds, where it laterally accelerates 
the outer stellar layers. This material, gravitationally conflned to the white dwarf, races 
around the star and, in ~2 s, converges at a point opposite to the location of the bubble's 
breakout, creating conditions capable of igniting the nuclear fuel and triggering a detonation 
that can incinerate the white dwarf and result in an energetic explosion. 

Of interest to the present study is the recent work of Kasen & Plewa (2005), who inves- 
tigate the spectral, and spectropolarimetric, consequences of the GCD model. They focus 
their investigation on the interaction of the expanding ejecta with an ellipsoidal, metal-rich 
extended atmosphere formed from the bubble of deflagration products (taken to be 57% 
Si, 27% S, 7.1% Fe, and 2.7% Ca, plus smaller amounts of other metals; see Khokhlov 
et al. 1993), and find that a dense, optically thick pancake of metal-rich material is formed 
at potentially large velocity on the side of the ejecta where the bubble emerged. For low 
atmosphere masses (e.g., resulting from a bubble of mass O.OOSMq), the pancake of mate- 
rial spans the velocity range 17,000-28,000 km s~^ and is geometrically detached from the 
bulk of the SN ejecta. This might explain the detached, high-velocity Ca II near-infrared 
(IR) triplet absorption seen in pre-maximum spectra of some SNe la (e.g., SN 200 lei, see 
§ 2.3.3). For larger bubble and, hence, atmosphere masses (matm ^ O.OI6M0), the absorb- 
ing pancake moves at lower velocities (e.g., 10,000-21,000 km s~^ for matm = O.O8M0) and 
could blend with the region of IMEs in the SN ejecta and potentially increase the blueshift 
of several of the spectral features. This could provide an orientation-dependent explana- 
tion for the origin of HV SNe la, whereby sight lines in which the pancake more completely 
blocks the photosphere produce the anomalously large velocities. The spectropolarimetric 
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consequences of this model are discussed in § 2.3.1. 

2.3. Spectropolarimetric Properties of SN la 

2.3.1. Supernova Polarization Mechanisms 

The first definitive proof that some SNe are polarized came from observations of the 
Type II SN 1987A (Cropper et al. 1988), which exhibited a modest temporal increase in 
continuum polarization during the photospheric phase, as well as sharp polarization mod- 
ulations across strong P-Cygni flux lines (Jeffery 1991b). While intrinsic polarization has 
now been estabhshed in over a dozen SNe (for recent reviews, see Wheeler 2000; Filippenko 
& Leonard 2004; Leonard & Fihppenko 2005), including at least three SNe la, the exact 
origin of both continuum and line polarization remains controversial. 

There are essentially two mechanisms by which supernova continuum polarization is 
thought to be produced: (1) a globally aspherical photosphere and electron-scattering atmo- 
sphere (e.g., Shapiro & Sutherland 1982; Hoflich 1991; Jeffery 1991a; Leonard et al. 2000a; 
Wang et al. 2001), and (2) ionization asymmetry produced by the decay of asymmetrically 

distributed radioactive ^^Ni, perhaps flung out into (Chugai 1992; Hoflich et al. 2001) or 
beyond (Kawabata et al. 2002; Leonard et al. 2002a) the expanding ejecta in clumps. The 
simplest, and most well-studied, globally aspherical geometry is that of an ellipsoid, and 
we shall make frequent reference to the "ellipsoidal model" in the following discussion. In 
the second model, in which an ionization asymmetry is present, continuum polarization 
is generated by light from the (either spherical or aspherical) photosphere scattering off 
of asymmetrically distributed free electrons that exist in clouds surrounding clumps of ra- 
dioactive ^^Ni. Variations on both of these polarization mechanisms are also possible. For 
instance, an aspherical distribution of ^^Ni could also result in SN polarization by providing 
an asymmetry in a source of optical photons (produced by the thermalization of 7-rays) 
relative to the scattering medium. 

To explain polarization modulations seen across spectral lines, it is important to differ- 
entiate between the emission peaks and blueshifted absorption troughs that are character- 
istic of P-Cygni proflles in total-flux spectra of SNe. Line peaks have usually been assumed 
to consist of intrinsically unpolarized photons. This is because although resonance scat- 
tering by a line is an inherently polarizing process (Jeffery 1991a), directional information 
for scattered photons is lost in an SN atmosphere since the timescale for randomizing coUi- 
sional redistribution of the relative level populations within the fine structure of the atomic 
levels of a line transition is much shorter than the characteristic timescale for absorption 
and reemission in a line (Hoflich et al. 1996). The assumption of intrinsically unpolarized 
emission lines is also commonly used to derive the interstellar polarization (ISP; § 2.3.2). 
Note, though, that if ionization asymmetry exists above the photosphere (from, e.g., clumps 
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of radioactive ^^Ni) , then even intrinsically unpolarized emission- line photons may become 
polarized by electron scattering within the SN atmosphere. 

For SNe la, line-blanketing, due largely to Fe, is particularly severe at wavelengths 
below ~5000 A, where theoretical models (Howell et al. 2001) suggest that nearly complete 
depolarization of any "continuum" light may be assumed. Conversely, the broad spectral 
region 6800 ^ A < 7800 A is largely devoid of line opacity (Kasen et al. 2004; Howell et al. 
2001), at least near maximum light, and thus may give a good indication of true continuum 
polarization level. 

The explanation of polarization changes in absorption troughs is controversial. One 
well-studied hypothesis, which follows logically from the ellipsoidal model, is that selective 
blocking of more forward-scattered and, hence, less polarized, light in P-Cygni absorption 
troughs results in trough polarization increases (Jeffery 1991a; Leonard et al. 2000b, 2001), 
potentially producing "inverted" P-Cygni polarization profiles when coupled with emission 
peak depolarization. Some SN modelers have claimed that this "geometrical dilution" 
mechanism, however, is rather poor at polarizing SN la light, and that polarization decreases 
should actually be seen in absorption troughs, since much of the light reaching the observer 
in those spectral regions has been absorbed and reemitted by the line (Howell et al. 2001). 

A key point is that, in either case, to produce line trough polarization changes under the 
ellipsoidal model, some continuum polarization must exist. In fact, Leonard h Filippenko 
(2001) show that the strength of a line trough polarization feature can be used to place a 
lower bound on the true intrinsic continuum polarization level under the ellipsoidal model 
according to 
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where Icont is the interpolated value of the continuum flux at the location of the line trough, 
Itrough is the total flux at the line's flux minimum, and Aptot is the total polarization change 
observed in the line trough. Thus, a critical test of the ellipsoidal model is to see whether the 
continuum polarization observed in the spectral region 6800-7800 A is suflicient to explain 
observed line trough polarization changes. Lack of signiflcant continuum polarization in 
an object with strong line features argues against the ellipsoidal model. An additional 
prediction is that no rotation of the polarization angle (PA) should be seen across the line, 
since the continuum and line-forming regions share the same geometry. 

A final way to produce a polarization change in a line trough is through asymmetry in 
the distribution of elements in the ejecta material located above the photosphere along the 
line-of-sight (1-o-s), hereafter referred to as the "clumpy ejecta" model. Asymmetry in the 
distribution of material with significant optical depth may unevenly block the underlying 
photospheric light, thereby producing a net polarization change and/or PA rotation in a line 
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trough, even when the photosphere is sphericaL The GCD model discussed in § 2.2 provides 
a natural mechanism by which to generate line-trough polarization through this model, 
since it predicts an optically thick pancake of high-velocity, metal-rich material overlying 
the photosphere on the side in which the bubble emerged. In particular, if the GCD 
model provides the correct mechanism by which to produce HV SNe la, then polarization 
changes in the absorption troughs of the strong, HV metal lines in these objects should be 
particularly pronounced since the 1-o-s necessarily intercepts a substantial fraction of the 
pancake. 

For SNe la, an interesting alternative to the simple ellipsoidal model for producing 
both line and continuum polarization is presented by Kasen et al. (2004), who explore the 
polarization consequences of a conical hole in the ejecta due to the interaction with a 
companion star, which we shall refer to as the "ejecta-hole" model. By considering various 
viewing angles and hole sizes, Kasen et al. (2004) demonstrate that both continuum and line 
polarization and PA changes can be generated. In general, viewing angles almost directly 
down the hole yield low continuum polarization with polarization increases in strong line 
troughs, whereas sight lines more nearly perpendicular to the hole result in larger continuum 
polarization (prominently seen in the spectral region 6800-7800 A) and "inverted P-Cygni" 
line polarization profiles; this latter case is nearly indistinguishable from the predictions of 
the simple ellipsoidal model. 

Ionization asymmetry is also capable of generating both continuum polarization as 
well as PA and polarization level changes through hne features, since both continuum and 
hne photons will scatter off of concentrations of free electrons in (or beyond) the ejecta. An 
impressive example of this mechanism potentially being at work is given by Chugai (1992) 
for the case of SN 1987A, in which the spectropolarimetry data, as well as asymmetries in 
the flux line profiles, are convincingly reproduced by the effects of two clumps of ^^Ni in 
the far (receding) hemisphere of the ejecta. 

2.3.2. Removing Interstellar Polarization 

A problem that plagues interpretation of all SN polarization measurements is proper 
removal of the ISP. Since directional extinction resulting from aspherical interstellar dust 
grains aligned by some mechanism along the 1-o-s to an SN can contribute a large polar- 
ization to the observed signal, an attempt must be made to remove it prior to analyzing 
SN spectropolarimetry data. This is notoriously difficult, although a number of different 
techniques have been advanced over the years, which we here summarize. 

An excellent way to derive Galactic ISP is through observations of distant, intrinsically 
unpolarized, "probe stars" close to (within ~0.5°) the 1-o-s to the SN (e.g., Leonard et al. 
2002a,b). Deriving the total ISP, which includes the contribution from dust in the host 
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galaxy, however, is more difficult. 

Tlie most basic tecfmique is to place limits on the ISP from reddening considerations. 
Since tlie same dust tliat polarizes starlight should redden it as well, it seems logical to 
expect a correlation between reddening and ISP. Because the alignment of dust grains is 
not total (or has multiple preferred orientations due to non-uniformity of the magnetic 
field along the 1-o-s), and grains are probably only moderately elongated particles, it is 
not surprising that, through the analysis of thousands of reddened Galactic stars, only 
an upper bound on the polarization efficiency of Galactic dust can be derived (Serkowski 
et al. 1975): ISF/Eb-v < 9.0% mag~^. However, the polarization efficiency of the dust 
in external galaxies is not well studied, and in one of the few investigations carried out 
to date, Leonard et al. (2002b) find compeUing evidence for polarization efficiency well in 
excess of the empirical Galactic limit for dust in NGC 3184 along the 1-o-s to the Type 
II-P SN 1999gi. It is not clear at this point whether meaningful constraints can thus be 
placed on the ISP of extragalactic SNe that are significantly reddened by host-galaxy dust. 
Nonetheless, total reddening arguments are still often used to set "reasonable" limits on 
the ISP. 

Another method to get a handle on the ISP is to assume axisymmetry for any SN 
asphericity, a situation that reveals itself through a straight-line distribution of points 
when the spectropolarimetry is plotted in the q-u plane. If axisymmetry exists, the ISP is 
constrained to lie along the axis defined by the line (Howell et al. 2001). Its absolute value, 
however, is uncertain without additional input. 

A more precise technique relies on the theoretical expectation that unblended emission 
lines consist of unpolarized light, and that any polarization observed in emission-line pho- 
tons comes from the ISP (Jeffery 1991b; Tran et al. 1997; Wang et al. 1996; Leonard et al. 
2001). This assumption is thought to be most vahd at early times, when any ^^Ni concen- 
trations are hkely to be below the photosphere. Note, though, that to use this method, 
care must be taken to isolate the emission-line photons from the underlying continuum 
hght (e.g., Tran et al. 1997), which may be intrinsically polarized. 

A related method is to assume that a particular spectral region is intrinsically com- 
pletely unpolarized, with all of the observed polarization therefore coming from ISP. Some 
empirical support exists that, in some objects, this may be the case for very strong emis- 
sion lines (Kawabata et al. 2002; Leonard et al. 2002a; Wang et al. 2004). For SNe la in 
particular, it has sometimes been assumed that the far blue spectral region (e.g., typically 
< 5000 A) satisfies this criterion due to the effect of the heavy line-blanketing and multiple, 
depolarizing, line scatters, largely due to iron-group elements. This qualitative expectation 
is demonstrated quantitatively by Howell et al. (2001), who present model polarization 
spectra resulting from delayed-detonation models in reahstic SN la atmospheres. More 
recently, Wang et al. (2005a) update this approach by choosing only specific spectral re- 
gions at blue wavelengths (A < 5000 A) that are not dominated by obvious individual fine 
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features in either flux or polarization for the ISP determination, rather than the arbitrary 
blue edge of the spectrum, as adopted by Howell et al. (2001). The central idea here is 
that while strong individual line features might impart their own polarization signature 
(e.g., through geometrical dilution in ellipsoidal models, or through blocking of the parts 
of the photosphere in clumpy-ejecta models), the spectral regions in between specific fea- 
tures probably decrease any effective continuum polarization by the numerous overlapping 
spectral lines at blue wavelengths. 

Finally, if polarimetry is obtained after the electron-scattering optical depth of the 
atmosphere has dropped well below unity, then it may be adequate to assume the observed 
polarization to be due entirely to ISP across the whole spectrum (e.g., Wang et al. 2003). 

Improperly removed, ISP can increase or decrease the derived intrinsic polarization, 
and it can change "valleys" into "peaks" (or vice versa) in the polarization spectrum. Since 
it is so difficult to be certain of accurate removal of ISP, it is generally safest to focus on (a) 
temporal changes in the polarization with multiple-epoch data, (b) distinct line features in 
spectropolarimetry having high signal-to-noise ratio (S/N), and (c) continuum polarization 
unlike the characteristic "Serkowski-law" wavelength dependence imparted by dust (e.g., 
Whittet et al. 1992). 

2.3.3. Previous SN la Polarimetry Studies 

Evidence for intrinsic polarization in SNe la was initially difficult to find, as the first 
investigations detected only marginally significant polarization among normal-brightness 
events observed near maximum, p < 0.2% (e.g., Wang et al. 1996, 1997). Broad-band 
polarimetry of one overluminous event observed over a month after maximum also found 
no intrinsic polarization down to a level of p ^ 0.3% (Wang et al. 1996). Significant 
advances have been made in the last few years. Leonard et al. (2000b) reported the first 
convincing, albeit weak, features in the polarization of an SN la, SN 1997dt, which was 
likely a subluminous event; these data are presented and analyzed in more detail in the 
present study. 

The first thorough spectropolarimetric study of an SN la is that by Howell ct al. (2001). 
In spectropolarimetry obtained at maximum light, the subluminous SN 1999by exhibits a 
polarization change of ~0.8% from 4800 A to 7100 A, and a sharp polarization modulation 
of ~0.4% across the strong Si II A6355 absorption. These features are explained within the 
context of an eUipsoidal model with a global asphericity of ~20%, observed equator-on. 
This physical picture was achieved by assuming the ISP to be the observed polarization 
at the blue edge of the spectrum, which was about 0.2% at 4800 A. With this choice of 
ISP, the inferred intrinsic polarization rises from 0% to about 0.8% from blue to red, with 
a sharp depolarization from 0.4% to near 0% across the Si II A6355 feature. An argument 
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supporting both the eUipsoidal model as well as the ISP choice is that after ISP removal 
the PA becomes nearly independent of wavelength. Howell et al. (2001) note, on the other 
hand, that a wavelength-independent PA also results if one assumes the ISP to be the 
observed polarization of the far red edge of the spectrum; in this case, the polarization 
modulation across the Si II A6355 becomes a polarization increase. Such a scenario could 
result from selective blocking of forward-scattered hght, as described in § 2.3.1. However, 
Howell et al. (2001) find the theoretical arguments supporting the former ISP to be more 
compelling. With this ISP, the redward rise in intrinsic polarization is attributed to the 
decreasing importance of line opacities, and the increased influence of continuum electron 
scattering at longer wavelengths. The depolarization across the Si II A6355 absorption is 
also attributed to the depolarizing effect of line scattering. Despite the detection of intrinsic 
polarization for SN 1997dt and SN 1999by, their unusual (subluminous) nature nevertheless 
left some doubt about intrinsic polarization in normal SNe la. 

That doubt has recently been put to rest with the work of Kasen et al. (2003) on 
SN 2001el. For this normal-luminosity event, the percent polarization changed from blue to 
red by ~0.4% in spectropolarimetry obtained 1 week before maximum brightness. However, 
the extraordinary feature here is the existence of distinct high- velocity Ca II near-infrared 
(IR) triplet absorption {v = 18,000-25,000 km s~^) in addition to the usual, lower-velocity 
Ca II feature. A similar, but much weaker, high-velocity feature had been previously 
observed in SN 1994D, and perhaps in other SNe la as well; the number of pre-maximum 
spectra covering the near-IR spectral range is small. The polarization is seen to change 
dramatically in this feature, by ~0.4%. The Ca II near-IR feature is examined by Kasen 
et al. (2003) in an elegant study, which concludes that it is likely due to photospheric 
obscuration by a clumped shell of high-velocity material. Using multi-epoch data, Wang 
et al. (2003) demonstrate that the nature of the polarization changes over the course of 
two weeks following this early epoch, becoming nearly undetectable a week after maximum 
brightness, further solidifying the case for intrinsic polarization at early times. Although not 
commented on by either study, it appears that the Si II A6355 line also shows a polarization 
modulation in the earliest epoch, and with a different PA from the Ca II near-IR feature. 
Although it does not affect the main results, a cautionary note on the difficulty of ISP 
determination is set by the fact that the two studies arrive at quite different values: Kasen 
et al. (2003) derive an ISP of very nearly 0% by assuming the blue edge of the spectrum 
to be unpolarized, whereas Wang et al. (2003) obtain an ISP of ~0.6% by attributing the 
observed polarization 38 days after maximum light, when SN 2001el is argued to be in the 
nebular phase, entirely to ISP. 

Most recently, Wang et al. (2005a) analyze a single epoch of pre-maximum spec- 
tropolarimetry of SN 2004dt, an HV SN la for which an ISP of qisp = 0.2 ± 0.1% and 
i^isp = —0.2 ± 0.1% is derived from the observed polarization of a handful of narrow, blue 
spectral regions in which no single spectral feature dominates in the total-flux spectrum. 
This results in rather low continuum polarization, p < 0.4%, but very strong modulations 



-13- 



across spectral lines. The polarization spikes reach 2% in the deep troughs due to Si II 
A4130 and A6355; lesser peaks are observed in features identified with Mg II A4471, a blend 
of Si II AA5041,5056 and Fe II AA4913, 5018, 5169, and the Ca II near-IR triplet. All line 
polarization has similar directional behavior in the q-u plane, suggesting a common origin. 
Interestingly, whereas other strong line features in the total-flux spectrum are characterized 
by strong polarization modulations, O I A7774 shows no polarization signature. These fea- 
tures are explained in terms of optically thick bubbles of IMEs, the result of partial burning, 
that are asymmetrically distributed within an essentially spherical oxygen substrate that 
remains from the progenitor. Note that Wang & Wheeler (1997) find polarization variation 
at a level of > 0.5% across the strong Si II A6355 feature of another HV SN la, SN 1997bp, 
in unpublished spectropolarimetry obtained near maximum light. 

Prom this small sample, the observations thus far suggest that normal and, perhaps, 
overluminous events are weakly polarized {p < 0.4%), with subluminous ones possessing 
somewhat larger values {p ^ 0.8%). HV SNe la have modest continuum polarization, 
but possess the highest line polarization of all, achieving polarizations of up to 2% in the 
strongest features. We assess the robustness of these tentative trends with the present 
study of four SNe la. 

2.4. The Type la Supernovae 1997dt, 2002bf, 2003du, and 2004dt 

SN 1997dt was discovered (Qiao et al. 1997) by the Beijing Astronomical Observatory 
Supernova Survey (Li et al. 1996) on 1997 November 22.44 (UT dates are used throughout 
this paper) at an unfiltered magnitude of ~15.3 in the Sbc galaxy NGC 7448. Images of 
the same field taken eight days earher show no star at the position of the SN to a limiting 
unfiltered magnitude of about 18.5. An optical spectrum obtained immediately (0.06 days) 
after discovery showed it to be a Type la event (Qiao et al. 1997); a subsequent examination 
by Li et al. (2001b) estimates the age at discovery to be — 7 ± 5 days relative to maximum 
light, which places the date of maximum at very roughly 1997 November 29. Tonry et al. 
(2003) report that SN 1997dt suffers from a host extinction of Ay = 0.46 mag, the median 
of the values derived from MLCS and "Baycsian Adapted Template Match Method" (J. 
L. Tonry et al. 2005, in preparation; Tonry et al. 2003) analyses of its unpublished light 
curves. An MLCS2k2 fit (S. Jha et al., in preparation) suggests that it is subluminous, 
with A = 0.94 (corresponding to a Ami5(i?) ~ 1.8 mag; see § 2.1), although there is a long 
shoulder of probability to lower A values, indicating that there is a wide range of light-curve 
shapes that can fit the sparse number of photometric points (2 in B, 3 in V, and 2 in 7). 
From a pre-maximum spectrum posted at the Center for Astrophysics' Recent Supernovae 
Page,^ we estimate TZ{Si II) 0.3 (see Nugent et al. 1995), which imphes Amir,{B) Ri 1.33 



^http://cfa- www.harvard.edu/cfa/oir/Research/supernova/RecentSN. html 
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mag from the correlation derived by Benetti et al. (2004), also consistent with a somewhat 
subluminous classification. 

SN 2002bf was discovered (Martin & Li 2002) by the Lick Observatory and Tenagra 
Observatory Supernova Searches (LOTOSS; Schwartz et al. 2000) on 2002 February 22.30 
at an unfiltcred magnitude of ~17, very close to the nucleus of the Sb galaxy PGC 029953. 
An image of the same field taken twenty days earlier showed nothing at the position of 
SN 2002bf to a limiting unfiltered magnitude of ~19. Optical spectra obtained on 2002 
March 6.21 by Mathcson et al. (2002) and on 2002 March 7.41 by Filippenko et al. (2002) 
identified it as a Type la event. Both groups noted that the expansion velocity derived from 
the absorption minimum of the Si II A6355 A line was significantly greater than normal for 
an SN la near maximum fight, indicating that it may be an HV SN la. A "spectral feature 
age" (Ricss ct al. 1997) of ± 2 days was derived from the March 6 spectrum (Matheson 
et al. 2002). The MLCS2k2 analysis of the light curves presented in § 4.2.1 yields 2002 
March 4.37 ± 0.50 as the date of maximum B light, consistent with the age on March 6 
derived by Matheson et al. (2002). 

SN 2003du was discovered (Schwartz & Holvorcem 2003) by LOTOSS on 2003 April 
22.4 at an unfiltered magnitude of ~15.9 in the SBd galaxy UGC 9391. An image of the 
same field taken fifteen days before discovery showed nothing at the position of SN 2003du 
to a limiting unfiltered magnitude of ~19. An optical spectrum obtained shortly thereafter, 
on 2003 April 24.06, identified it as a Type la event roughly two weeks before maximum 
fight (Kotak et al. 2003). The MLCS2k2 analysis of the fight curves presented in § 4.2.1 
yields 2003 May 6.12 ± 0.50 as the date of maximum B light, consistent with the earfier 
spectral age estimate. It also agrees with the epochs of maximum determined by the recent 
photometric studies of SN 2003du by Anupama et al. (2005) and Gerardy et al. (2004). 

SN 2004dt was discovered (Moore & Li 2004) by the Lick Observatory Supernova 
Search (Filippenko et al. 2003) at an unfiltered magnitude of ~16.1 in the SBa galaxy 
NGC 799 on 2004 August 11.48. An image of the same field taken ten days earlier showed 
nothing at the position of SN 2004dt to a limiting unfiltered magnitude of ~18. Spectra 
taken within 2 days of discovery by Gal- Yam (2004), Patat et al. (2004), and Salvo et al. 
(2004) confirmed it to be an SN la before maximum light. Patat et al. (2004) noted 
that several absorption lines showed high expansion velocities, a result confirmed by the 
recent study by Wang et al. (2005a), suggesting that, like SN 2002bf, SN 2004dt is an 
HV SN la. A preliminary analysis of the light curves of SN 2004dt shows that maximum 
fight occurred near 2004 August 20 (W. Li, personal communication). A series of Hubble 
Space Telescope (HST) UV spectral observations was obtained as part of program GO-10182 
(P.I. Filippenko), and will be analyzed in a future paper. 
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3. Observations and Reductions 

We obtained single-epoch spectropolarimetry of SN 1997dt, SN 2002bf, SN 2003du, and 
SN 2004dt on days 21, 3, 18, and 4 (respectively) after maximum light. We also obtained 
additional optical spectroscopy and BVRI photometry of SN 2002bf and SN 2003du. For 
SN 2002bf, our photometry samples —10 to 57 days from the time of maximum hght, with 
one additional flux spectrum taken on day 9 after maximum. For SN 2003du, our ground- 
based photometry covers —5 to 113 days from maximum, with one additional epoch on day 
436 taken using the High Resolution Channel (HRC) of the Advanced Camera for Surveys 
(ACS) on board HST. Five additional spectral epochs sample its development from days 
24 to 82 after maximum. 



3.1. Photometry 

3.1.1. Ground-Based Photometry of SN 2002bf and SN 2003du 

All ground-based photometric data were obtained using either the 0.76-m Katzman 
Automatic Imaging Telescope (KAIT; Fihppenko et al. 2001; Li et al. 2003) or the Nickel 1 
m reflector (Li et al. 2001a), both located at Lick Observatory. Figures 1 and 2 show KAIT 
images of PGC 029953 and UGC 9391, the host galaxies of SN 2002bf and SN 2003du, 
respectively. Also labeled in the KAIT images are the "local standards" in both fields that 
were used to measure the relative SN brightness on non-photometric nights. We obtained 
12 epochs of Johnson- Cousins BVRI photometry (Johnson et al. 1966, for BV\ Cousins 
1981, for RI) for SN 2002bf, all taken with KAIT, and 38 epochs of BVRI photometry 
for SN 2003du, 33 of them taken with KAIT and five with the Nickel telescope. We also 
obtained three pre-maximum unfiltered observations of SN 2002bf with KAIT (approxi- 
mating the it! band, see Li et al. 2003), and one additional epoch of RI photometry with 
KAIT of SN 2003du. 

For the photometry we employed the usual techniques of galaxy "template" subtraction 
(Li et al. 2000, and references therein), point-spread function fitting (Stetson 1991, and 
references therein), and using "local standards" to determine the BVRI brightnesses of the 
SNe on non-photometric nights; in general, we closely followed the technique detailed by 
Leonard et al. (2002d). We note that the galaxy subtraction procedure for SN 2002bf was 
particularly challenging since it is only 4'.'1 from its host galaxy's center. 

The absolute cahbration of the SN 2002bf field was accomphshed on the photometric 
nights of 2002 May 14 and 2004 March 17 with the Nickel telescope, and 2003 February 
3 and 2004 March 18 with KAIT, by observing several fields of Landolt (1992) standards 
over a range of airmasses in addition to the SN 2002bf field. The absolute calibration of the 
SN 2003du field was similarly derived from data taken on the photometric nights of 2003 
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May 31, June 1, June 26, and August 27 with the Nickel telescope, and of 2003 May 22 and 
2004 March 18 with KAIT. The color terms used to transform the filtered instrumental 
magnitudes to the standard Johnson-Cousins system are those of Foley et al. (2003). We 
list the measured BVRI magnitudes and the la uncertainties, taken as the quadrature sum 
of a typical photometric error and the la scatter of the photometric measurements from 
all of the photometric nights, of the local standard stars in Tables 1 and 2. 

After deriving the BVRI magnitudes of the SNe based on a comparison with each 
of the local standards, we took the weighted mean of the individual estimates as the final 
standard magnitude of the SNe at each epoch in each filter. The results of our ground-based 
photometric observations are given in Tables 3 and 4 and shown in Figures 3 and 4. The 
reported uncertainties come from the quadratic sum of the photometric errors (reported 
by DAOPHOT) and the transformation errors. For SN 2002bf, the uncertainty produced 
by the difficult galaxy-subtraction process contributed the majority of the error on nights 
with low S/N. 

3.1.2. Hubble Space Telescope Photometry of SN 2003du 

We obtained HST images during the course of two orbits of a 29" x 26" field of 
view centered on SN 2003du on 2004 July 15, 434 days after B^^^,, with the ACS/HRC 
detector through filters F435W, F555W, F625W, and F814W (hereafter referred to as 
B,V,R, and/, respectively), as part of our Snapshot survey program (GO-10272; P.I. 
Li) to investigate the late-time photometric behavior and environment of nearby SNe. 
SN 2003du was detected in all images. Total exposure times in BVRI were, respec- 
tively, 1680 s (data archive designation j8z441011/3011), 960 s (j8z442011/4011), 720 s 
(j8z441021/3021), and 1440 s (j8z442021/4021). 

To derive the i/ST photometry, we followed as closely as possible the procedure detailed 
by Sirianni et al. (2005), including correction for the effects of SN light contaminating the 
background region, aperture corrections, the "red-halo" effect (for the /-band), and CTE 
degradation (Riess 2004). We translated the resulting instrumental magnitudes to the 
standard Johnson- Cousins BVRI system by using the coefficients and color corrections 
tabulated by Sirianni et al. (2005). The final results of our //ST photometry are included 
in Table 4 and displayed in Figure 4. 

3.2. Spectropolarimetry and Spectroscopy 

We obtained single epochs of spectropolarimetry for SN 2002bf and SN 2003du on 2002 
March 7 and 2003 May 24, respectively, with the Low- Resolution Imaging Spectrometer 
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(Oke et al. 1995) in polarimetry mode (LRISp; Cohen 1996)^ at the Cassegrain focus of 
the Keck-I 10-m telescope. We observed SN 1997dt on 1997 December 20 with LRISp 
using the Keck II 10-m telescope, and SN 2004dt on 2004 August 24 with the Kast double 
spectrograph (Miller & Stone 1993) with polarimeter at the Cassegrain focus of the Shane 
3-m telescope at Lick Observatory. We reduced the polarimetry data according to the 
methods outlined by Miller et al. (1988) and detailed by Leonard et al. (2001) and Leonard 
& Filippenko (2001). 

The polarization angle offset between the half- wave plate and the sky coordinate system 
was determined by observing the following polarized standard stars from the list of Schmidt 
et al. (1992) and setting the observed F-band polarization position angle (i.e., 9v, the 
debiased, flux-weighted average of the polarization angle over the wavelength range 5050- 
5950 A; see Leonard et al. 2001) equal to the cataloged value: BD +64°106 (1997 December 
20), BD +59°389 (2002 March 7), and HD 161056 (2003 May 24). On the night of 2004 
August 24, we averaged the polarization angle offsets derived from observations of three 
polarized standards from the Schmidt et al. (1992) list, HD 204827, BD +59°389, and 
HD 19820; the individual offsets were internally consistent to within 1°. To check for 
instrumental polarization, the following null standards taken from the lists of Turnshek 
et al. (1990), Mathewson & Ford (1970), Schmidt et al. (1992), and Berdyugin et al. 
(1995), were also observed: HD 94851 (1997 December 20), HD 57702 (2002 March 7), 
HD 109055 and BD +32°3739 (2003 May 24), and HD 212311 (2004 August 24). All stars 
were measured to be null to within 0.1%, which is also our estimate of the systematic 
uncertainty of a continuum polarization measurement made with either the Keck or Lick 
spectropolarimeters (e.g., Leonard et al. 2001). 

Additional specifics of the observations of SN 2002bf taken on 2002 March 7, including 
an investigation of the potential impact of second-order light contamination and instru- 
mental polarization (both shown to be minimal) in the setup used on this night, are given 
by Leonard et al. (2002a). Second-order light contamination is not a concern for our obser- 
vation of SN 1997dt due to its limited spectral range. For SN 2003du, the use of a dichroic 
(D560) to split the beam near 5600 A eliminates second-order light contamination on the 
red side. Our spectropolarimetric observation of SN 2004dt was taken in a setting that 
included the use of an order-blocking filter (GG455) to prevent contamination by second- 
order fight at red wavelengths. Beyond ~9000 A, however, second-order contamination 
may exist, but for reasons similar to those discussed by Leonard et al. (2002a) we believe 
it to have minimal impact for this particular object. 

To derive the total-fiux spectra, we extracted all one-dimensional sky-subtracted spec- 
tra optimally (Horne 1986) in the usual manner. Each spectrum was then wavelength 
and fiux calibrated, and was corrected for continuum atmospheric extinction and telluric 



^Instrument manual available at littp://www2. keck.hawaii.edu/inst/lris/pol_quickref.html. 
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absorption bands (Wade & Home 1988; Bcssell 1999; Matheson et al. 2000). With the 
exception of the spectropolarimetric observations of SN 1997dt and SN 2002bf, all spectra 
were taken near the parallactic angle (Filippenko 1982), so the spectral shape should be 
quite accurate. Table 5 lists the spectropolarimetric and spectral observations for all four 
SNe. Figures 5-8 show the observed spectropolarimetry data of the four objects, and Fig- 
ures 9 and 10 show the complete series of spectra obtained for SN 2002bf and SN 2003du, 
respectively. 

4. Analysis 

4.1. Spectroscopy 

Our spectrum of SN 1997dt, taken ~21 days after maximum, shows typical features for 
an SN la at this phase (Fig. 5a). Similarly, our spectral sequence of SN 2003du (Fig. 10) 
closely follows the evolution of the normal-luminosity SN la 1994D (Patat et al. 1996; 
Filippenko 1997) in terms of the strengths and blueshifts of line features. This is consistent 
with the analyses of Anupama et al. (2005) and Gerardy et al. (2004), in which spectra taken 
before and shortly after maximum light were also examined. This convincingly establishes 
SN 2003du as a spectroscopically "typical" SN la (Branch et al. 1993). 

As discussed in § 2.4, the spectra of both SN 2002bf and SN 2004dt are peculiar in 
one regard: the blueshifts of many of the spectral lines, most noticeably Si II A6355, occur 
at significantly higher velocity than is typical for an SN la at this phase, and indicate that 
these are both HV SNe la. 

For comparison, we have measured the velocities of the Si II A6355 line in a number of 
other HV SNe la from our database, along with the spectroscopically normal SN 1994D, the 
subluminous SN 1991bg (e.g., Filippenko et al. 1992a), and the overluminous SN 1991T 
(e.g., Fihppenko et al. 1992b). We present the results in Table 6 and Figure 11, from 
which it is clear that both SN 2002bf and SN 2004dt belong to the class of HV SNe la; 
indeed, SN 2002bf is the most extreme HV SN la yet observed for its epochs. The figure 
also suggests that line velocity is not strongly correlated with luminosity, although the 
expansion velocity of the subluminous SN 1991bg is somewhat lower than typical values. 

Figure 12 presents a spectral comparison near maximum light of three HV SNe la 
(SN 2002bf, SN 2002bo, and SN 2004dt) with the spectroscopically normal SN 1994D. The 
excessive blueshift of the Si II A6355 trough is obvious for the three HV SN la compared 
with SN 1994D. While it must be cautioned that the spectra span a range of ages of about 
five days near maximum light, a time when significant spectral development occurs, the 
blueshift differences are much greater than can be explained by age differences alone. In 
addition, the Si II line is significantly stronger in the HV SNe la compared with SN 1994D, 
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with equivalent widths of > 140 A compared with ~100 A for SN 1994D. Conversely, the 
O I A7774 absorption is relatively weaker in the HV SN la events, with equivalent widths 
of < 100 A in the HV SNe la compared with 125 A for SN 1994D. The fact that the Si II 
line is significantly stronger in the HV SNe la sample, and the O I line relatively weaker, 
is consistent with the scenario in which a greater fraction of C and O is burned to IMEs 
in HV SN la than in more typical events; it also follows naturally from the GCD model, 
since the obscuring pancake is formed from the products of oxygen burning. 

While HV SNe la share many spectral characteristics, it is clear that they also exhibit 
spectral diversity. For instance, whereas the absorption trough of the Ca II near-IR triplet is 
significantly blueshifted for SN 2002bf and SN 2002bo relative to SN 1994D (16, 300 km s~^ 
and 14, 500 km s~^ for SN 2002bf and SN 2002bo, respectively, compared with ~10,900 km 
s~^ for SN 1994D, where we have assumed Aq = 8579 for the Ca II near-IR triplet, a value 
derived using the prescription given by Leonard et al. 2002c), it has a more normal blueshift 
{v = 10,900 km s~^) in SN 2004dt. The equivalent widths of the Ca II near-IR lines also 
show a suggestive trend, with both SN 2002bf and SN 2002bo having widths of > 200 A, 
while both SN 2004dt and SN 1994D have equivalent widths of -100 A. This may indicate 
that the explosive nucleosynthesis did not proceed up to Ca as far out in the atmosphere 
(or in the bubble in the GCD scenario) of SN 2004dt as it did in the other two HV SNe. 

Finally, it is clear from Figure 12 that not all lines share the extreme velocities seen 
in the Si II A6355 feature. For instance, as first noticed by Wang et al. (2005a) in a pre- 
maximum spectrum of SN 2004dt, the S II "W" feature in the spectra of SN 2002bo and 
SN 2004dt indicates velocities comparable to those in SN 1994D; for SN 2002bf they are 
somewhat larger, but still not near the velocity of the Si II A6355 fine. Wang et al. (2005a) 
propose that this may indicate that sulfur is more confined to the lower-velocity, inner 
region. 

4.2. Photometry 

4-2.1. Ground-Based Photometry 

The results of the MLCS2k2 apphcation to the photometry of SN 2002bf and SN 2003du 
are given in Table 7. For details of the MLCS2k2 procedure used, see S. Jha et al. (in 
preparation); an overview of the technique is provided by Riess et al. (2005). The MLCS2k2 
analysis finds SN 2002bf to be of typical luminosity. SN 2003du is slightly overluminous, 
although its pre-maximum spectral evolution (Anupama et al. 2005) demonstrates that it 
is not a SN 1991T-like event, as the strength of the Si II A6355 feature is comparable to 
that seen in spectra of normal SNe la. 

For SN 2002bf, E^ax occurred on 2002 March 4.37 ± 0.50, which is 10 days after 
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discovery and 5 days before our filtered observations commenced. For SN 2003du, -Bmax 
occurred on 2003 May 6.i2 ± 0.50, wliicli is 14 days after discovery, and the same day that 
our BVRI observations began. Our derived date of maximum hght for SN 2003du agrees 
with those found by Anupama et al. (2005) and Gerardy et al. (2004) using independent 
data sets. 



4.2.2. HST Photometry of SN 2003du 

Late-time photometry (e.g., t > 200 d) of SNe la exists for only a handful of objects 
(see, e.g., Milne et al. 2001, and references therein). From the small sample, there are two 
main features of note. First, SN la decline rates are typically much faster than the decay 
slope of ^^Co — > ^^Fe of 0.98 mag (100 d)~^ predicts. This decay mechanism is thought to 
be primarily responsible for powering the luminosity from the early nebular phase out to 
~1000 days for SNe of all types. Essentially, the ^^Co decays release most of their energy 
in the form of 7-rays, which, given enough optical depth, can become trapped in the ejecta 
and Compton scatter off free electrons. The energetic electrons generate optical photons 
primarily through ionization and excitation of atoms, and the ejecta are transparent to 
these photons. A small portion (~3.5%, see Arnett 1979) of the total ^^Co decay energy 
comes in the form of positrons, which may deposit their kinetic energy in the ejecta and 
then annihilate with electrons, producing two 7-ray photons of energy = 171^0^. The 
steeper decline seen in the late-time photometry of SNe la has been explained by significant 
transparency of the ejecta to 7-ray photons (Milne et al. 1999) and positron escape (Milne 
et al. 2001). 

A quantitative measure of the late-time decline is given by Cappellaro et al. (1997), who 
investigate how the decline from maximum to 300 days after peak V brightness, denoted 
^?7T.3oo(l^), correlates with intrinsic SN brightness as derived from the Ami5{B) parameter. 
Although limited by small sample size (5 objects), they find a convincing correlation, with 
Am3oo(^) going from 6.7 mag to 8.4 mag as the sample runs from overluminous (SN 1991T) 
to subluminous (SN 1991bg) events, with the normal-brightness SN 1994D characterized 
by Am3oo(V^) = 7.3 mag. A second feature that has thus far been seen in only two SN la 
events (SN 1991T and SN 1998bu; see Schmidt et al. 1994 and Cappellaro et al. 2001, 
respectively) is a sudden flattening of the late-time optical light curves, which has been 
attributed to the contribution of a light echo from foreground dust clouds. 

Our photometric data from HST, taken 436 days after i?max, allow us to investigate the 
late-time photometric behavior of SN 2003du. From the inset of Figure 4, it is clear that 
SN 2003du, hke all SNe la analyzed before it, declines significantly faster than the ^^Co — > 
^^Fe decay rate. In fact, we measure the average decay rate in V (shown by Milne et al. 2001 
to track the bolometric luminosity of an SN la quite accurately) of SN 2003du between our 
last two photometric epochs on days 113 and 436 to be AV — 1.47 ± 0.02 mag (100 d)~^. 
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The slope also appears to have been rather constant throughout the period between our 
two last epochs, as indicated by the good agreement of the late-time data taken from 
Anupama et al. (2005) near day 300 and the decay slope determined from our data alone. 
Using the Anupama et al. (2005) data point and our estimate of V"max (Table 7), we derive 
^rnsoo^V) = 6.74, which is most consistent with the values found by Cappellaro et al. 
(1997) for overluminous events. That SN 2003du may be somewhat overluminous was also 
suggested by its peak V magnitude of — 19.67 ± 0.02 mag, which is 0.17 mag brighter than 
the fiducial template used in the MLCS2k2 procedure (Table 2). There is no evidence from 
our data of any contribution to the SN brightness from a light echo, although we note that 
the major indications of additional contributions to the apparent brightness of SNe 1991T 
and 1998bu did not become obvious until epochs > 500 d. 



4.3. Reddening 

4.3.1. Techniques to Estimate SN la Reddening 

Accurate determination of SN reddening is crucial both for deriving intrinsic SN prop- 
erties as well as interpreting spectropolarimetry, since the same dust that reddens SN light 
can also polarize it as discussed in § 2.3.2. When multi-band photometry is available, the 
MLCS2k2 technique can accurately estimate the total extinction (see Riess et al. 2005 for 
discussion). When this is lacking, other methods must be used. 

Galactic extinction along the 1-o-s is accurately estimated by the dust maps of Schlegel 
et al. (1998, hereafter SFD) to an estimated precision of ~15%. For host-galaxy extinction, 
a rather crude approach, to which many investigators resort in the absence of photometry, 
is to employ the rough correlation found between the total equivalent width (Wx) of the 
interstellar (IS) Na I D doublet (AA5890, 5896) and reddening (Barbon et al. 1990). The 
Barbon et al. correlation has been subsequently improved upon by Munari & Zwitter 
(1997), who derive a more precise relation that uses just the equivalent width of the Na I 
D2 (A5890) line. Both relations, however, warrant healthy degrees of skepticism since 
sodium is known to be only a fair tracer of the hydrogen gas column (especially in dense 
environments, where sodium may be heavily depleted; e.g., Cohen 1973), from which the 
dust column is then estimated. The dust-to-gas ratio also varies significantly among galaxies 
(e.g., Issa et al. 1990), and it could well be that in the case of SNe, circumstellar, rather than 
interstellar, dust is present, for which the dust-to-gas ratio (or the extinguishing properties 
of the dust itself) could be unusual. A final complication is that, at the resolution typical of 
most optical spectra, individual absorption components along the 1-o-s, whose contributions 
should be considered separately to determine the total reddening, are blended into a single 
profile. In such situations, the Munari & Zwitter (1997) relation formally yields an upper 
limit to the reddening, and caution must be used since the derived value could seriously 
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overestimate the actual reddening. Nonetheless, the Munari & Zwitter (1997) relation is 
stiU frequently used to get approximate reddening values, or upper limits, especially in the 
case of null detections of IS Na I D lines. 



4.3.2. The Reddening of the Four SNe la 

For SN 2002bf and SN 2003du we shall adopt the MLCS2k2 values of E{B-V) = 
0.08 ± 0.04 mag and E{B—V) — 0.01 ± 0.01 mag (respectively) reported in Table 7 (for 
Ry = 3.1). For SN 1997dt, Tonry et al. (2003) report a host-galaxy reddening value 
of E{B-V) Host = 0.15 mag (assuming Ry = 3.1), the median of the values derived from 
MLCS2k2 and "Bayesian Adapted Template Match Method" analyses. The SFD dust maps 
predict E{B-V)f^ = 0.06 mag for SN 1997dt, for a total estimate of E{B-V) = 0.21 
mag. 

Our spectrum of SN 1997dt also exhibits strong Na I D IS absorption at both the 
redshift of the host galaxy and the MW. The resolution of our spectrum, ~5 A (Table 5), 
while sufficient to deblend the Na I D doublet, is not fine enough to resolve the individual 
absorption components that likely contribute to the Dl and D2 profiles. For host-galaxy 
absorption, we measure Wl°^ = 0.77 A, with W^^ = 0.42 A and W^^ = 0.35 A. This 
translates to E(B-V)^% < 0.21 mag from the Munari & Zwitter (1997) relation. For 
the MW, we measure 1^^°* = 0.76 A, with W^^ = 0.44 A and W^^ = 0.32 A, yielding 
E{B—V)^^j^ < 0.23 mag from the Munari & Zwitter relation. The values given by SFD and 
Tonry et al. (2003) are consistent with the upper limits derived from the sodium relation. 
The purported accuracy of the SFD Galactic reddening value and the upper limit set by the 
Munari & Zwitter relation for the host reddening would indicate an upper reddening limit 
of E{B—V)tota.i < 0.27 mag. This limit will prove to have important implications when 
we examine the spectropolarimetry of SN 1997dt in § 4.4.2. The fact that the Munari 
& Zwitter relation overpredicts the reddening for both the host galaxy and, especially, 
the MW, may indicate that multiple, unresolved components make up the IS line profiles. 
Given the upper reddening limit, we conclude E{B—V) — 0.21 ± 0.06 mag for SN 1997dt. 

For SN 2004dt, E{B-V)f^ = 0.03 mag, and Wang et al. (2005a) report a range of 
E{B—V)tota.i — 0.14 to 0.2 mag from analysis of the color of their pre-maximum spectrum. 
However, they also note that SN 2004dt does not show noticeable Na I D IS lines at 
the redshift of NGC 799. We confirm the lack of Na I D lines in our spectrum as well, 
although its poor resolution (~18 A, see Table 5) makes deriving even an upper limit to 
the strength of the Na I D IS lines difficult. Formally, our procedure yields W\(3a) = 0.1 A 
for host-galaxy Na I D, which translates to a predicted upper limit of E{B—V)^H\^ < 0.02 
mag. We thus confirm the discrepancy noted by Wang et al. (2005a) between the lack of 
IS sodium absorption and reddening inferred from other methods. It could be that the 
spectroscopic peculiarities of HV SNe la produce photometric irregularities that affect the 
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general reddening relations. As mentioned in § 2.2, Benetti et al. (2004) do find photometric 
peculiarities, albeit minor, in their study of another HV SN la, SN 2002bo. However, since 
the Na I D relation, especially for the host galaxy, is also prone to error, we have no 
convincing way to decide between the two estimates. We thus take the simple average of 
the low and high reddening estimates, and incorporate the disparity into our final estimate's 
uncertainty. This yields E{B—V) — 0.11 ± 0.06 mag as our best reddening estimate for 
SN 2004dt. 



4.4. SpectropolEirimetry of Four SNe la 

To summarize: SN 1997dt is likely a somewhat subluminous event that is thought to 
be reddened by E{B-V) = 0.21 ± 0.06 mag, with an upper hmit of 0.27 mag. SN 2002bf 
and SN 2004dt are HV SNe la, with SN 2002bf being the most extreme example yet 
observed. Finally, SN 2003du is a shghtly overluminous event that is minimally reddened. 
At maximum light, a previous subluminous event has been found to be moderately polarized 
{p ^ 0.8%), whereas normal to overluminous examples have been less so {p < 0.4%). The 
two HV SNe la that have reported polarization measurements show the greatest polarization 
features yet observed for any SN la type, with values approaching 2% in the strongest lines 
of SN 2004dt (§ 2.3.3). 

We now examine our single-epoch spectropolarimetry of SN 1997dt, SN 2002bf, SN 
2003du, and SN 2004dt obtained on days 21, 3, 18, and 4, respectively, after maximum light. 
For each SN, we shall first discuss the observed polarization and then attempt to remove 
the ISP through the technique of Wang et al. (2005a), which assumes that spectral regions 
lacking strong individual flux or polarization features at blue wavelengths (A < 5000 A) 
are intrinsically unpolarized. Since the choice of the specific spectral regions to use for ISP 
determination is admittedly somewhat subjective, we shall be careful to point out how our 
conclusions would change with other ISP choices. 

Following ISP removal and examination of the data in the q-u plane, we calculate the 
intrinsic polarization of the SN by rotating the q-u axes through a single angle, 9, that 
places the greatest degree of polarization change across the spectrum along the rotated q 
axis. The angle by which the axes are rotated is determined through a uniform- weight, 
least-squares fit to the ISP-subtracted data in the q-u plane. The polarization degree 
measured along the rotated q axis is then referred to as the rotated Stokes parameter 
(RSP; Tran 1995; Leonard et al. 2001), whereas that measured along the rotated u axis is 
denoted URSP.^ 



^Note that RSP and URSP, as defined here, are identical to the Wang et al. (2001) definitions of the 
polarization along the "dominant" and "orthogonal" axes, Pd and Pg, respectively. That is, by definition, 
RSP shows polarization strength variation along the chosen polarization angle (i.e., the dominant axis), 
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We shall then analyze the resulting spectropolarimetry within the context of the models 
described in § 2.3.1. By construction, the greatest degree of polarization change, especially 
from the continuum regions, should occur in the RSP. Examination of URSP is especially 
important in the line features, however, as it can be used to discriminate between the 
predictions of the ellipsoidal and clumpy-ejecta models. In particular, the simple ellipsoidal 
model demands the existence of a single value of 9 capable of placing all line polarization 
changes along the RSP (i.e., URSP should show no polarization features). 

44.1. Two HV SNe la: SN 2002hf and SN 2004dt 

We begin by considering the spectropolarimetry of the two HV SNe la, SN 2002bf and 
SN 2004dt, which are displayed in Figures 6 and 8, respectively. Both objects show very low 
levels of observed polarization, with pv = 0.03%, 9v = 62° for SN 2002bf and pv = 0.25%, 
9v — 146° for SN 2004dt, where pv and 9v approximate the rest-frame F-band polarization 
and polarization angle derived by calculating the debiased, flux-weighted averages of q and u 
over the interval 5050-5950 A (see Leonard et al. 2001). The polarization of SN 2004dt may 
exhibit an upward trend with wavelength, increasing from about 0.3% at blue wavelengths 
to ~1.0% at the red edge (A = 9600 A); a similar but smaller trend exists in the data for 
SN 2002bf. Neither object shows significantly different polarization in the "continuum" 
region 6800-7800 A from that observed in the heavily line-blanketed region below 5000 A 
(§ 2.3.1), suggesting that the intrinsic continuum polarization is quite low for both objects. 

The salient features of both data sets are the extraordinarily large polarization mod- 
ulations (~2%) across certain P-Cygni lines, most notably the Ca II near-IR triplet for 
SN 2002bf and the Si II A6355 line for SN 2004dt. Smaller features may be discerned 
across other fines in both objects. 

Before attempting ISP removal, it is instructive to compare our SN 2004dt data with 
those obtained by Wang et al. (2005a) taken eleven days earlier, wfien tfie SN was about 
seven days before maximum light. The interval separating the two observations is one 
marked by rapid spectral and photometric evolution. For SN 2004dt, the velocity of the 
minimum of the high-velocity Si II A6355 fine recedes by about 3700 km (from 17, 200 



whereas the URSP shows polarization strength variation along the axis that is orthogonal to the dominant 

axis in the q-u plane. The main difference between our construction of the RSP in the observed polar- 
ization plots (Figures 5-8) and the ISP-subtracted data (Figures 15, 17, and 19) is that in the former the 
polarization angle about which the RSP is determined is a smoothly varying function of wavelength (i.e., 
it is the polarization angle, 0, smoothed over many bins; see Leonard et al. (2001) for more details on the 
procedure), whereas in the latter, RSP and URSP are determined with respect to a single, imchanging 
polarization angle. For consistency with spectropolarimetry work in other fields, as well as our own prior 
studies, we shall continue to use the RSP and URSP designations, rather than the Pd and Po designations 
of Wang et al. (2001). 



- 25 - 



to 13,500 km s~^), and the weaker, lower-velocity lines such as S II AA5612,5654 (Aq = 
5635 A) decrease by about 2100 km s^^ (from 11,000 to 8900 km s^^). Remarkably, in 
spite of the great changes that occur in the flux spectrum between the two epochs, we 
find the spectropolarimetry to be virtually unchanged. The overall polarization, in both 
magnitude and polarization angle, is consistent with the earlier epoch, and the detected 
line-polarization features, most notably those at Si II A6355 and the Ca II near-IR triplet, 
are also very similar. As with the Wang et al. data, our observations do not show significant 
polarization modulation across the important O I A7774 line, although this line is now much 
weaker in the total-flux spectrum than it was earlier. Weaker polarization features, which 
are clearly detected in the earlier epoch, are difficult to confirm in our lower S/N data, but 
our data are not inconsistent with the earlier measurements, within the errors. 

As discussed in § 2.3.3, Wang et al. (2005a) explain the spectropolarimetry of SN 2004dt 
in terms of optically thick bubbles of IMEs that are asymmetrically distributed within an 
essentially spherical oxygen substrate that remains from the progenitor material. Wang 
et al. (2005a) further note that the polarization behavior of the high-velocity lines is sim- 
ilar to the low-velocity lines in the early epoch, implying that the structures that obscure 
the photosphere have great radial extent. Within this context, then, it may not be surpris- 
ing that the polarization features have not evolved significantly between the two epochs, 
despite the great evolution of hne velocity. 

In order to directly compare the spectropolarimetry of SN 2002bf and SN 2004dt 
with each other, we must attempt to remove the ISP. To estabhsh the ISP, we apply the 
technique of Wang et al. (2005a) and choose the spectral regions indicated in Figures 6a 
and 8a. For SN 2002bf, this yields (gisp,iiisp) = (0.01%, 0.05%), or ISPi„£« = 0.05% at 
Oisp = 39° at an assumed peak wavelength of the ISP of Amax = 5500 A. For SN 2004dt we 
derive (q'isp, Misp) = (0.3%, —0.2%), which is within the uncertainty of the ISP found by 
Wang et al. (2005a) through this same approach, {qisp,Uisp) = (0.2%, —0.2%). In order to 
facilitate direct comparisons between our data and those of Wang et al. (2005a), we shall 
adopt their ISP value for our study of SN 2004dt as well. 

After removal of the small amount of ISP for SN 2002bf, we derive a best-fitting axis 
with PA = 123° (Fig. 13), about which we calculate the intrinsic RSP and URSP shown 
in Figure 15. In a similar way, for SN 2004dt we derive a best-fitting axis of PA = 146°, 
which may be compared with the PA of ~150° found by Wang et al. (2005a). Again, 
since the two values do not differ significantly, we adopt the Wang et al. direction for ease 
of comparison between the two data sets. The axis and ISP choices for SN 2004dt are 
indicated in Figure 14, and the resulting RSP and URSP are shown in Figure 15. 

From examination of Figure 15, it is clear that the two events have spectropolarimetric 
similarities. Both show low overall polarization and modulations across the Si II A6355 
and Ca II near-IR triplet absorptions. The detailed character of the hne polarizations do 
differ, however. For SN 2002bf there is a 2% polarization change across the Ca II feature. 
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with a more modest (~0.4%) feature detected in the Si II line. (When discussing overall 
"polarization change" across a line feature without specific reference to either the RSP or 
URSP, we effectively mean the quadrature sum of the changes seen in the two parameters.) 
For SN 2004dt, the situation is reversed. In the case of the Ca II line, this may be related 
to its relative strength and velocity in the two spectra, as both the equivalent width and 
velocity of the feature in the flux spectrum are much greater in SN 2002bf than in SN 2004dt. 
As discussed earlier (§ 4.1), it is possible that burning to Ca occurred more extensively in 
SN 2002bf than it did in SN 2004dt. The explanation for the Si II line disparity, however, 
is not so obvious, as the lines have similar strengths and velocities. 

Clearly, strong line and weak continuum polarization levels disfavor the simple ellip- 
soidal model. Applying Eq. (1) to the data for SN 2002bf and SN 2004dt yields lower 
bounds on the expected continuum polarization of Pcont ^ 0.7% and Pcont > 1-3%, re- 
spectively. For both objects, however, we measure p < 0.4% in both the observed and 
ISP-subtracted data for the spectral region 6800-7800 A, which is largely devoid of line 
opacity in SNe la near maximum (§ 2.3.1). The ellipsoidal model therefore appears to be 
ruled out as the explanation for the polarization characteristics of these objects. 

The clumpy-ejecta model, on the other hand, provides a natural explanation for high 
line and low continuum polarization levels (§ 2.3.1). In particular, the GCD scenario can 
successfully explain (1) the high line velocities, (2) the large polarization change in the line 
troughs, (3) the great radial extent of the obscuring material (i.e., for matm = O.O8M0, the 
obscuring pancake spans the range 10,000-21,000 km s~^), and (4) the lack of significant 
continuum polarization. A remaining challenge is to explain the lack of any polariza- 
tion change across the O I A7774 line, especially in the earfier Wang et al. (2005a) data for 
SN 2004dt, when this line is extremely strong in the total-fiux spectrum. Recent nucleosyn- 
thesis calculations based on multi-dimensional (2D and 3D) hydrodynamical simulations 
of the thermonuclear burning phase in SNe la show that as much as 40% to 50% of the 
ejected matter in SNe la is unburned carbon and oxygen (Travaglio et al. 2004). This has 
led Wang et al. (2005a) to propose that the primordial oxygen is nearly spherically dis- 
tributed, within which asymmetrically distributed IME clumps, or an absorbing pancake 
in the GCD model, are embedded. 

With this in mind, an important prediction of the clumpy-ejecta model (including the 
GCD scenario) is that the O I A7774 line should, in fact, show a polarization change with a 
polarization angle that differs by 90° from what is observed in the Si II and Ca II lines since 
its distribution is essentially the inverse of these IMEs. However, some oxygen is probably 
also contained in the clumps or pancake, since it is produced by explosive carbon burning 
in small quantities and should be present when such burning products as magnesium exist, 
whose spectral signature is unequivocally seen in the spectra. This would tend to reduce 
the polarization level in the O I A7774 line. Neither our data, nor those of Wang et al. 
(2005a), are of sufficiently high S/N to detect such a change, and it must be left to future. 
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higher-S/N studies focused especially at early times when the O I A7774 feature is strong. 
Finding such a PA change would further strengthen the case for the clumpy-ejecta and/or 
GCD scenario. 

We note that the cjccta-holc model of Kasen et al. (2004) is also capable of producing 
large line polarization with weak continuum polarization for sight-lines near to the hole 
(§ 2.3.1). Arguing against this model in these cases, though, is that it offers no natural 
explanation for why high line velocities should be associated with sight lines near to the 
hole. In fact, the Kasen et al. models predict lower absorption velocities when viewing 
down the hole. 

In all, then, our study finds that HV SNc la arc, as a group, characterized by much 
stronger line-polarization features than are seen in other SN la varieties. The ellipsoidal 
model is incapable of explaining the polarization characteristics of these objects, whereas 
the clumpy-ejecta and ejecta-hole models are more successful. On balance, the case for 
clumpy ejecta appears to be the most convincing explanation, with the GCD model in- 
vestigated by Kasen & Plewa (2005) able to reproduce many of the observed spectral and 
spectropolarimetric features. 

44.2. SN 1997dt 

We next turn to the likely subluminous SN 1997dt, which was observed about three 
weeks past maximum light. Figure 5 reveals an extraordinarily high level of observed 
polarization, pv — 3.46% at 9v — 112°. This is by far the largest polarization yet observed 
for an SN la. A distinct polarization feature is detected in the Fe II A4555 trough and 
probably also in the Si II A5972 + Na I D and Si II A6355 lines. The degree of change in 
the Fe II A4555 feature reaches nearly 1% in the q parameter. 

This amount of observed continuum polarization is surprising, given that our best total 
reddening estimate predicts an upper bound on the ISP of only 1.89% from the Serkowski 
et al. (1975) relation (see § 2.3.2). If we trust both the reddening estimate and the upper 
ISP bound, then an intrinsic SN polarization of at least 1.57% must exist to explain the 
observed polarization, far higher than has been indicated for any previous SN la. However, 
the technique of Wang et al. (2005a) suggests a much higher ISP level: {qisp,uisp) = 
(2.53%, 2.68%), or ISP^nax = 3.60%, 9isp = 113°, for X^^^ = 6500 A, the wavelength of 
maximum ISP that yields the most convincing Serkowski-law fit to the observed polarization 
(Fig. 5d). If this truly is the ISP, then it impfies an extraordinarily high polarization 
efficiency for the dust along the 1-o-s to SN 1997dt: ISP/Eb-v ~ 18% mag~^, which is 
double the observed Galactic limit of 9% mag~^. Of course, some of this discrepancy could 
be removed if the true reddening were greater than we suspect. However, even allowing the 
reddening to equal the upper limit of E{B—V)totai < 0.27 mag set in § 4.3.2 still requires 
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the dust-polarization efficiency to exceed tlie Galactic limit. 

It thus appears that we face a stark choice: either SN 1997dt has the highest intrinsic 
polarization of any SN la yet observed, or the dust along the 1-o-s has an exceptionally 
high polarization efficiency. The epoch of our observation of SN 1997dt is unique for a 
subluminous SN la, so there is no empirical database from which to draw expectations and 
help decide between these two options. 

The situation for SN 1997dt, while perplexing, in fact is not unique: a similar palette 
of possibilities presented themselves in our previous study of a single epoch of spectropo- 
larimetry of SN 1999gi, an SN II-P that also had a low reddening and large observed 
polarization (Leonard & Filippenko 2001; Leonard et al. 2002d). Like SNe la, SNc II-P as 
a group have historically shown very low intrinsic continuum polarization. After consider- 
ing many polarization production mechanisms, including polarization due to newly formed 
dust in the SN ejecta and dust reflection by one or more off-center dust blobs external to 
the SN, we concluded that the most hkely explanation for the polarization of SN 1999gi 
was that the host-galaxy dust along the 1-o-s possesses a very high polarization efficiency, 
1SP/E[B—V] = 31^9^% mag~^, which remains the largest value yet inferred for a single 
sight line in either the MW or an external galaxy (Leonard et al. 2002b). 

There are arguments favoring a similar explanation here. First, a Serkowski law rea- 
sonably fits the observed continuum polarization of SN 1997dt (Fig. 5d). Further, at 21 
days past maximum, SN 1997dt may be nearing the end of its photosphcric phase, a time 
when spectropolarimetry may be losing its efficacy as an asymmetry probe due to a lack of 
electrons available to scatter the light. If we believe this to be the case, then wc should not 
expect large polarization, even if the SN is highly aspherical. If wc accept the large ISP, 
then we naturally would like to know whether it is due to dust in the MW or NGC 7448. 
Unfortunately, we have not observed any distant Galactic "probe" stars (§ 2.3.2) near to 
the 1-o-s to estimate the Galactic ISP. There are, however, reasons to suspect that it is low. 
First, 20 stars within 10° of the 1-o-s are listed in the agglomeration of stellar polarization 
catalogs by Heiles (2000), and the greatest observed polarization is only 0.3%. Second, the 
great majority of the reddening is due to host-galaxy dust, since E{B—V)jiost = 0.15 mag 
while E{B—V)mw = 0.06 mag (§ 4.3.2). As was the case with SN 1999gi, we would again 
conclude that it is the dust within the host galaxy that must have the extraordinarily high 
polarization efficiency. 

There are arguments to oppose this, however. Previous SN la polarization studies 
have found intrinsic polarizations increasing toward red wavelengths (§ 2.3.3), which could 
certainly mimic a Serkowski law over the limited wavelength band covered by our obser- 
vations. Further, the line polarization features demonstrate that the SN must possess at 
least some intrinsic polarization. Finally, the relatively late phase of the observation, in- 
voked previously to argue against high intrinsic polarization, can also be used to argue in 
favor of it: at the stage immediately before an SN begins the transition to the nebular 



-29- 



phase, the deepest layers of the ejecta are revealed. Although in need of confirmation by 
detailed modeling, for this epoch one can plausibly argue that an optical photosphere still 
exists with sufficient optical depth to electron scattering, perhaps even reaching the single- 
scattering limit, which is the most polarizing atmosphere (Hoflich 1991). If the explosion 
mechanism itself is asymmetric, the largest imprint of the asymmetry would presumably be 
in the innermost layers, which could lead to a very large intrinsic polarization that reveals 
itself just at this late phase. 

In fact, such an effect is seen in the spectropolarimctry of the SN II-P 2004dj as it 
transitions to the nebular phase (D. C. Leonard et al., in preparation). However, in the case 
of SN 2004dj, the strong increase is observed primarily in the "continuum" region 6800- 
7800 A (and in a few strong fine troughs) , not in the overall level across the whole spectrum. 
Although in need of confirmation by detailed modeling, depolarizing line blanketing may 
also be significant for SNe la at blue wavelengths at this epoch, which would again argue 
for a large ISP as the explanation of the high observed polarization. One possibility that 
circumvents this difficulty is that asymmetrically distributed concentrations of radioactive 
Ni, recently uncovered in the thinning ejecta at this relatively late epoch, are responsible 
for the large polarization. 

Curiously, the single, high-polarization observation of SN 1999gi occurred at a similar 
stage of its evolution, right at the end of the optical plateau that characterizes the photo- 
spheric phase in SNe II-P. It is thus unfortunate that no other spectropolarimetric epochs 
were obtained for either event (SN 1999gi and SN 1997dt) to serve as a basis for compari- 
son. Multi-epoch data covering the transition from the photospheric to the nebular phases 
of SNe of all types will certainly help reveal more definitively the physical explanation for 
such high observed polarizations at these late times. 

It may be possible to gain insight into the cause of the observed polarization from 
examination of the spectropolarimetry data in the q-u plane, shown in Figure 16. ISP 
originating from a single source (i.e., characterized by a single PA) will spread intrinsically 
unpolarized data points along a line in a direction that intersects the origin in the q-u 
plane. Finding that data lie predominantly along such a line, and exhibit a Scrkowski law 
spectral shape, provides compelling evidence that a large, single source of ISP dominates 
the observed signal. Such was the case for SN II-P 1999gi (Leonard et al. 2002b). Similarly, 
for SN 1997dt, an elongation of the data points from blue to red wavelengths in a direction 
that roughly points back toward the origin also exists. 

When the ISP derived earlier (ISP„iax = 3.60%, O^sp = 113°, for A^ax = 6500 A) 
is removed, the wavelength dependence of the polarization largely disappears (Fig. 16), 
strengthening the argument for a single, dominant ISP source. The data are then seen to 
lie along an axis with a fairly well-defined PA of 58° (Fig. 16), against which we calculate the 
RSP and URSP shown in Figure 17. Although of individually low to moderate significance, 
the polarization features in the Fe II A4555 and Si II A6355 troughs, seen in both RSP 
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and URSP, do seem to prefer the same general direction, with most of the modulation 
occurring along the URSP axis (e.g., the direction perpendicular to the main axis in the 
q-u plane). Taken at face value, consistent line-trough polarization changes in a direction 
differing from that preferred by the continuum favors an origin in the selective blocking of 
photospheric light by clumpy and asymmetrically distributed intermediate and iron-peak 
elements overlying the photosphere as opposed to an ellipsoidal scenario. 

Other plausible ISP choices, however, yield different conclusions. For instance, a 
smaller ISP of 2.6% at the same PA yields a more nearly spherical constellation of points 
centered near (g, u) = (—0.7%, —0.07%), with the hne excursions now pointing back toward 
the origin. Since some modelers predict that, in SNe la with an ellipsoidal asphericity, po- 
larization decreases may exist in absorption troughs (§ 2.3.1), the ellipsoidal model cannot 
therefore be ruled out in this the cause of the polarization. Given the marginal 

significance of the weaker features, additional interpretation of the line-trough polarization 
degrees and directions is probably not warranted with these data. 

In conclusion, we find evidence for a large ISP contribution to the observed polariza- 
tion of SN 1997dt, probably in the range 2.6% < ISP < 3.6%. This implies a polarization 
efficiency for the dust along the 1-o-s in NGC 7448 that exceeds the Galactic limit. We 
do, however, also find evidence for polarization intrinsic to the object, most convincingly 
in specific line features and, perhaps, in the continuum as well. A combination of asym- 
metrically distributed radioactive Ni and synthesized IMEs overlying the photosphere may 
provide the simplest explanation for the line and potential continuum polarization at this 
late phase, although the eUipsoidal model cannot be definitively ruled out. 

44.3. SN 2003du 

Our spcctropolarimctry of SN 2003du represents the highest S/N data of our study, and 
permits a more detailed analysis of line features than has been possible for the other objects. 
SN 2003du presents a low level of observed polarization across most of the spectrum, 
with pv — 0.04%, 9v — 17° (Fig. 7). The polarization exhibits an increasing trend with 
wavelength, rising from nearly zero at blue wavelengths to ~0.2% at the red edge of the 
spectrum. There are distinct and significant polarization changes across several absorption 
features, including the Ca II near-IR triplet, the Si II A6355 line, probably a few weaker 
lines such as Fe II A4924, and very tentatively the Ca II H & K absorption. 

Applying the technique of Wang et al. (2005a), and choosing the regions indicated in 
Figure 7a to estimate the ISP, yields (gisp,nisp) = (-0.02%, 0.0%), or ISP^na^ = 0.02%, 
^isp = 90°, for an assumed Amax = 5500 A. This very low ISP is consistent with the 
negligible reddening found earlier. Furthermore, on the same night the SN 2003du data were 
taken, we observed the distant Galactic star BD -|-50°1593 (spectral type F8, V = 10.64 
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mag), located just 0.28° from the 1-o-s of SN 2003du, and found it to be null to within 
0.1%. From its spectroscopic parallax, we estimate BD +50°1593 to be at least 190 pc 
away which, at a Galactic latitude of 53°, satisfies the criterion of Tran (1995) that a good 
MW "probe" star be more than 150 pc from the Galactic plane. We thus have multiple 
reasons to suspect httle ISP contaminating the data. 

After removal of the minimal ISP, a well-defined axis with PA = 107° is derived 
(Fig. 18), about which we calculate the intrinsic RSP and URSP, shown in Figure 19. 
Compared with what was seen in the HV SNe la, and even SN 1997dt, the line-polarization 
features are not large, amounting to no more than 0.3%. However, the very high S/N of 
these data makes the detections unequivocal, and establishes intrinsic polarization in a 
spectroscopically and photometrically normal SN la at the latest phase yet observed. 

The behavior of the line-trough polarization for the Ca II near-IR triplet and Si II 
A6355 are very similar. Both show sharp increases of ~0.2% in RSP, as well as overall 
increases in URSP. Modest RSP increases of < 0.1% may also be discerned in the Fe II 
A4924 and Ca II H & K absorptions. At this phase, the O I A7774 feature in SNe la is 
quite weak and, coupled with telluric A-band contamination, makes definitive detection 
of polarization modulation in this important region difficult; the observed changes are at 
about the level of the statistical noise. 

The commonality between the polarization behavior of the Si and Ca lines argues for 
similar origins. Given the high S/N of the data, we can examine rather closely the basic 
predictions of the ellipsoidal model that (a) the polarization angle should be independent of 
wavelength, and (b) the overall polarization should increase from blue to red wavelengths, 
with an expectation of p ^ 0% at A < 5000 A (if the line opacity remains strong in these 
regions at the epoch of observation; § 2.3.1). 

With our initial choice of ISP, the first criterion is clearly not satisfied, as there are 
obvious excursions in URSP throughout the spectrum, with especially sharp modulations 
occurring across the strongest lines. In fact, for ISP values constrained to lie along the 
symmetry axis we are unable to find any value that convincingly satisfies both criteria of 
the eUipsoidal model: values in the upper-right quadrant (e.g., [gisPj'Wisp] = [0.3%, 0.18%]) 
produce a polarization decrease with wavelength across the spectrum, while ISP values 
in the lower-left quadrant that satisfy the criterion of p — > 0% at blue wavelengths (e.g., 
[gisp, Misp] = [—0.1%, —0.07%]) result in strong PA changes across the spectral lines. Choos- 
ing an ISP point far beyond the constellation of data points, which may serve to make the 
PA changes less objectionable (although still statistically significant), has the unfortunate 
consequence of straining the limits implied from the very low reddening, at least for dust 
with normal polarizing efficiency. 

In addition, examining the URSP behavior of the Si II A6355 and Ca II near-IR lines 
in Figure 19 more closely, we see that the generally increasing trends in URSP across the 
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lines show sharp decreases right at the locations of peak RSP modulation, although the 
statistical significance of the modulations, especially for Ca II, is not high. If the abrupt 
changes in URSP are real, such structure is readily explained under the clumpy-ejecta model 
by variations in the distribution of the IMEs as a function of radius in the expanding ejecta. 
Chugai (1992) also demonstrates that such sharp changes in fine features can be produced 
by excitation asymmetry. On the other hand, there is no obvious mechanism to produce 
such an effect in the simple eUipsoidal models. We thus conclude that it is difficult to 
reconcile the basic predictions of the ellipsoidal model with the data for SN 2003du. 

We therefore suspect either clumps in the ejecta overlying the photosphere or ionization 
asymmetry as the cause of the inferred line and, perhaps, continuum polarization, and are 
led again to disfavor the ellipsoidal model as the cause of the intrinsic polarization of this 
SN. 



5. Conclusions 

We present post-maximum single-epoch spectropolarimetry of four SNe la, bringing 
to six the number of SNe la thus far examined in detail with spectropolarimetry during 
the early phases. The four objects span a range of spectral and photometric properties, yet 
all are demonstrated to be intrinsically polarized. This suggests that asphericity and/or 
asymmetry may be a ubiquitous characteristic of SNe la in the first weeks after maximum 
light. The nature and degree of the polarization varies considerably within the sample, 
but in a way that is consistent with, and extends, previously suspected trends. Our main 
spectropolarimetry results are as follows: 

1. SN 2002bf and SN 2004dt, both HV SNe la observed shortly after maximum bright- 
ness, exhibit the largest polarization features yet seen definitively for any subtype of 
SN la, with modulations of up to ~2% in the troughs of the strongest fines. The 
overall polarization level of both objects is minimal, at least at blue wavelengths; 
there is a possible trend of increasing polarization with wavelength for both objects, 
though neither shows significant polarization in the "continuum" region 6800-7800 A. 
The ISP contamination is not thought to be large in either object. 

2. SN 1997dt, believed to be a somewhat subluminous event, has the highest observed 
overall polarization of any SN la yet studied, pv = 3.46%, at 21 days past maximum 
light. This demands either an extraordinarily large polarization efficiency for the 
dust along the 1-o-s in NGC 7448, the largest intrinsic SN la polarization thus far 
found, or perhaps some combination of the two. The observed polarization rises by 
about 0.5% from blue (A = 4300 A) to red (A = 6700 A) wavelengths, approximating 
a Serkowski-law ISP curve rather convincingly, albeit one with a somewhat unusual 
peak wavelength (Amax ~ 6500 A). A polarization modulation of nearly 1% in the 
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strong Fe II A4555 absorption, and a more modest change of ~ 0.3% in the Si II A6355 
hne, demonstrate that the SN does possess intrinsic polarization features. However, 
we conclude that ISP is responsible for the bulk of the overall polarization that is 
observed, with 2.6% < ISP < 3.6%, and that the polarization efficiency of the dust 
along the 1-o-s in NGC 7448 likely exceeds the empirical Galactic limit. 

3. SN 2003du is a slightly overluminous SN la. Our spectropolarimetry, taken 18 days 
after maximum light, is the highest S/N data obtained for our sample of objects. 
It reveals a low continuum polarization that increases by ~0.3% from blue to red 
wavelengths, with distinct changes of ~0.2% detected in the Si II A6355 and Ca II 
near-IR triplet lines; smaller changes may be detected in weaker lines. The very 
similar behavior of the polarization in the two strongest lines, in both magnitude and 
direction in the g-i* plane, suggests a common polarization origin. ISP is thought to 
be minimal. 

Ordered by increasing strength of line-polarization features in SNe la, we find as fol- 
lows: ordinary/overluminous < subluminous < HV SNe la, with the strength of the line- 
polarization features increasing from 0.2% in the slightly overluminous SN 2003du to 2% 
in both HV SNe la in our study. Absolute continuum polarization levels are more difficult 
to establish, due largely to uncertainties in the ISP, but there arc compelling reasons to 
believe that at least three of our objects possess very little intrinsic polarization in spectral 
regions outside of specific fine features. The Howell et al. (2001) study of SN 1999by and 
our data on SN 1997dt provide some evidence that continuum polarization may be higher 
in subluminous objects than in other types. 

There are a number of alternatives from which to choose for the origin of polarization 
in SNe la, including global asphericity (e.g., the ellipsoidal model), ionization asymmetry, 
and clumps in the ejecta overlying the photosphere. The small, redward rise in the overall 
polarization level that is discerned in at least three of our objects can be reproduced 
by models possessing either global asphericity or an ionization asymmetry. Under the 
ellipsoidal model, the levels of continuum polarization for SN 2002bf, SN 2003du, and 
SN 2004dt imply minor to major axis ratios of around 0.9 if viewed equator-on (Hoflich 
1991; Wang ct al. 2003). This level of asphericity would produce a luminosity dispersion of 
about 0.1 mag for random viewing orientations (Hoflich 1991), which could explain some 
of the dispersion seen in the brightness-decline relation of SNe la. If the proposed global 
asphericity is more complicated, then the luminosity of a single SN la may depend on 
viewing angle in a non-trivial way such that, even for a large sample of objects, an overall 
bias to slightly higher or lower values may result. 

The potential for significantly greater continuum polarization, perhaps of in the 

likely subluminous SN 1997dt observed ~21 days after maximum would imply a more severe 
distortion, of at least 20%, from the models of Hofiich (1991). Since most cosmological 
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applications of SNe la rely on data acquired closer to maximum light, such a late-time 
asphericity, even if common, would not seriously affect the utility of SNe la as distance 
indicators. 

To explain the ubiquitous line polarization, the simple ellipsoidal model is effectively 
ruled out for three of our objects, including, most convincingly, the two HV events. From 
a number of lines of reasoning, the most convincing explanation is partial obscuration of 
the photosphere by clumps of newly synthesized IMEs forged in the explosion. 

For the HV SNe la, in particular, the GCD model studied by Kascn & Plewa (2005) 
provides a plausible explanation for many of the observed spectral and spcctropolarimetric 
characteristics. It predicts the existence of an optically thick pancake of material with sig- 
nificant radial extent that partially obscures the optical photosphere, producing larger line 
velocities and equivalent widths for many spectral features, and stronger line-trough polar- 
ization than is seen in more typical events. These qualitative expectations are borne out 
by our data: The line features of SN 2002bf and SN 2004dt possess the strongest polariza- 
tion modulations and greatest equivalent widths of our sample. The astonishing similarity 
between our epoch of spectropolarimetry of SN 2004dt, ~4 days after maximum, and that 
presented by Wang et al. (2005a) from 11 days earlier, provides compelling evidence that 
the obscuring material also possesses great radial extent in the thinning ejecta. 

That SNe la may be separable into different groups based on their spectropolarimetric 
characteristics yields one more clue to assist in narrowing down progenitor possibilities 
and/or models for the physics of the explosion. The assertion that some fraction of the IMEs 
in the ejecta of SNe la may be confined to bubbles or filaments is, however, a rather blunt 
discriminatory tool: At the present stage of theoretical modeling, dcfiagration, delayed- 
detonation, off-center dclayed-detonation, and GCD models can all plausibly be argued to 
produce clumps in the ejecta (e.g., Wang et al. 2005a, and references therein). The specific 
predictions of the GCD model in particular need to be further examined, preferably in 
full three-dimensional simulations, to test whether the explosion mechanism itself is viable, 
and, if so, whether it can quantitatively reproduce the observed characteristics of at least 
some SNe la in detail. On the observational front, higher S/N data, preferably obtained at 
multiple epochs, will help to narrow down the possibilities as well. With the steady advances 
being made in the theoretical understanding of these events, and the growing rate of SNe la 
studied in detail with spectropolarimetry, prospects for improving our understanding of 
these events are bright. 
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Fig. 1.— fi-band image of PGC 029953 taken on 2002 March 15 with KAIT, with the local 
standards listed in Table 1 marked. For clarity the inner 0'5 of PGC 029953 is expanded 
with the contrast adjusted to clearly show the SN (inset). SN 2002bf (SN) is measured to 
be 0'/5 east and 4'.'1 north of the low-surface-brightness center of PGC 029953. 



-42- 




Fig. 2.— V-hand image of UGC 9391 taken on 2003 May 19 with KAIT, with the local 
standards listed in Table 2 marked. SN 2003du (SN) is measured to be 9'.'3 west and 14'.'9 
south of the center of UGC 9391. 
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Fig. 3. — BVRI light curves for SN 2002bf from Table 3. For clarity, the magnitude scales 
for BRI have been shifted by the amounts indicated. The open symbols on days —10, —8, 
and —4 represent KAIT unfiltered magnitudes translated to the R band, as described by 
Li et al. (2003). Arrow indicates the epoch of our spectropolarimetric observation. 
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Fig. 4. — BVRI light curves for SN 2003du from Table 4. For clarity, the magnitude scales 
for BRI have been shifted by the amounts indicated. In most cases the error bar is smaller 
than the plotted symbol. Arrow indicates the epoch of our spectropolarimetric observation. 
The inset shows the complete light curve, including the HST ACS/HRC observations on 
day 436. Dotted lines connect the day 113 and day 436 data points, and open symbols are 
data from Anupama et al. (2005) on day 301 (2004 January 3). 
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4500 5000 5500 6000 6500 
Rest Wavelength (A) 

Fig. 5. — Observed polarization data for SN 1997dt obtained on 1997 Dec. 20, when the 
SN was about 21 days past maximum B brightness, (a) Total flux, in units of 10^"*^^ ergs 
cm~^ with prominent absorption features identified. Line identifications in this 

and all figures are from Li et al. (2001a), and references therein. Horizontal dashes indicate 
spectral regions used to determine the ISP. (b-c) Normalized q and u Stokes parameters. 
(d) Observed degree of polarization, (e) Polarization angle in the plane of the sky. The 
total-flux spectrum is shown at 2 A bin~^, whereas the polarization data are binned at 20 A 
bin~^ to improve the S/N. The displayed polarization p is the "rotated Stokes parameter" 
in this and all figures; see discussion in § 4. The thin lines in b, c, and d show a Serkowski 
law ISP characterized by the parameters derived in § 4.4.2. Note that the NASA/IPAC 
Extragalactic Database (NED) recession velocity of 2194 km s~^ for NGC 7448 has been 
removed in this and all figures showing spectra of SN 1997dt. 
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Fig. 6. — As in Figure 5, except for SN 2002bf obtained 2002 Mar. 7, when the SN was 3 
days past maximum B brightness. The total-flux spectrum is shown at 2 A bin^^, whereas 
the polarization data are binned at 100 A bin^-^ to improve the S/N. Note than many of the 
large jumps seen in 9 are an artifact of the data going through the origin in the q-u plane. 
The thin lines in b, c, and d show a Serkowski law ISP characterized by the parameters 
derived in § 4.4.2. The NED recession velocity of 7254 km s"^ for PGC 029953 has been 
removed in this and all figures displaying spectra of SN 2002bf. 
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Fig. 7.— As in Figure 5, except for SN 2003du obtained 2003 May 24, about 18 days after 
Bmax- The total-flux spectrum is shown at 2 A bin~^, whereas the polarization data are 
binned at 50 A bin~^ to improve the S/N. Note that the region 5400-5700 A was adversely 
affected by the dichroic used to split the beam, and thus could not be calibrated with 
certainty. The thin lines in b, c, and d show a Serkowski law ISP characterized by the 
parameters derived in § 4.4.2. The NED recession velocity of 1914 km s"^ for UGC 9391 
has been removed in this and all figures displaying spectra of SN 2003du. 
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Fig. 8. — As in Figure 5, except for SN 2004dt obtained 2004 Aug. 24, about 4 days after 
Bmax- The total-flux spectrum is shown at 5 A bin^^, whereas the polarization data are 
binned at 50 A bin^^ to improve the S/N. The thin lines in b, c, and d show a Serkowski 
law ISP characterized by the parameters derived in § 4.4.2. The NED recession velocity of 
5915 km s~^ for NGC 799 has been removed in this and all figures displaying spectra of 
SN 2004dt. 
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Fig. 9. — Optical flux spectra of SN 2002bf labeled by epoch with respect to B maximum. 
For clarity, the spectra have been scaled and shifted vertically by arbitrary amounts. Note 
that while the spectrum at t = 3 days is binned 2 A bin~^, the spectrum at 7 days has 
been binned at 20 A bin~^ in order to increase the S/N. 
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Fig. 10. — Optical flux spectra of SN 2003du labeled by epoch with respect to B maximum. 
For clarity, the spectra have been scaled and shifted vertically by arbitrary amounts. 
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Fig. 11. — Blueshifts of Si II A6355 for nine SNe la, from Table 6. Lines connecting data 
points for each SN are shown to help guide the eye, and are not the results of formal fits. 
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Fig. 12.— Spectra of four SNe la (from top down: SN 2002bf, SN 2002bo, SN 2004dt, 
SN 1994D) with time from -Bmax indicated, demonstrating the diversity of P-Cygni 
absorption-hne blueshifts. All but SN 1994D are considered to be "high-velocity" SNe la, 
due to their unusually broad and highly blueshifted Si II A6355 A line. The expected lo- 
cations of several additional line features, if blueshifted by 15,400 km s~^ (the blueshift of 
the Si II A6355 line for SN 2002bf) from their rest wavelengths, are indicated by vertical 
lines. The spectra have been dereddened by the following amounts: E{B—V) = 0.08 
mag (SN 2002bf; § 4.3.2), E{B-V) = 0.43 mag (SN 2002bo; Benetti et al. 2004), 
E{B-V) = 0.11 mag (SN 2004dt; § 4.3.2), and E{B-V) = 0.06 mag (SN 1994D; Patat 
et al. 1996). 
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Fig. 13. — Observed polarization data for SN 2002bf from Figure 6 plotted in the Stokes 
q-u plane. Each point represents a bin 100 A wide, with the symbol size increasing from 
the blue end of the spectrum (A ~ 3900 A) to the red (A 8500 A). The error bars 
show the typical statistical uncertainty of a data point. The large, open triangle near the 
origin indicates the adopted ISP, and the thick solid line indicates the axis about which 
the intrinsic RSP and URSP were calculated following ISP removal, with the arrowhead 
defining the (arbitrarily chosen) direction of increasing RSP. 
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Fig. 14. — Polarization data for SN 2004dt from Figure 8 plotted in the Stokes q-u plane. 
Each point represents a bin 50 A wide, with the symbol size increasing from the blue end 
of the spectrum (A ^ 4400 A) to the red (A 9600 A). The error bars show the typical 
statistical uncertainty of a data point. The large, open triangle indicates the adopted 
ISP, and the thick solid line indicates the axis about which the intrinsic RSP and URSP 
were calculated following ISP removal, with the arrowhead defining the (arbitrarily chosen) 
direction of increasing RSP. 
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Fig. 15.— Polarization data for the HV SNe la SN 2002bf and SN 2004dt following ISP 
removal and calculation of the RSP and URSP with respect to the axes indicated in Fig- 
ures 13 and 14, respectively. The thick, smooth lines without error bars in the RSP plots 
are the total-flux spectra, arbitrarily scaled and shifted for comparison of features. 
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Fig. 16. — Polarization data for SN 1997dt from Figure 5 plotted in the Stokes q-u plane, 
where filled circles represent the observed polarization, and filled squares represent the 
polarization data after removal of the ISP derived in § 4.4.2. Each point represents a bin 
20 A wide, with the symbol size increasing from the blue end of the spectrum (A 4300 A) 
to the red (A f» 6750 A). The two open triangles that are connected by a dotted line indicate 
the range of the ISP for our adopted ISPmax of 3.60% at ^ = 113° for Amax = 6500 A over 
the observed wavelength range; the smaller triangle is the ISP at the blue edge of the 
spectrum and the larger triangle is the ISP at Amax, which is very near the red edge of 
our spectrum. The error bars show the typical statistical uncertainty of a data point. The 
thick solid line indicates the axis about which the intrinsic RSP and URSP were calculated 
following ISP removal. The (arbitrarily chosen) direction of increasing RSP is indicated by 
the arrowhead. 
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Fig. 17. — Polarization data for SN 1997dt following ISP removal and calculation of the 
RSP and URSP with respect to the axis indicated in Figure 16. The thick smooth line in 
the RSP plot is the total-flux spectrum, arbitrarily scaled and shifted for comparison of 
features. 
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Fig. 18. — Polarization data for SN 2003du from Figure 7 plotted in the Stokes q-u plane. 
Each point represents a bin 50 A wide, with the symbol size increasing from the blue end 
of the spectrum (A ~ 3500 A) to the red (A 9300 A). The error bars show the typical 
statistical uncertainty of a data point. The opeM triangle indicates the adopted ISP, and the 
thick solid line indicates the axis about which the intrinsic RSP and URSP were calculated 
following ISP removal, with the arrowhead defining the arbitrary direction of increasing 
RSP. 
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Fig. 19. — Polarization data for SN 2003du following ISP removal and calculation of the 
RSP and URSP with respect to the axis shown in Figure 18. The thick smooth line in 
the RSP plot is the total-flux spectrum, arbitrarily scaled and shifted for comparison of 
features. 
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Table 1. Magnitudes of Local Standards: SN 2002bf 



Star Bias) V{av) R{aR) I{ai) 



1 12.640 (0.014) 12.226 (0.011) 11.986 (0.020) 11.709 (0.003) 

2 17.112 (0.021) 16.485 (0.010) 16.116 (0.006) 15.776 (0.027) 

3 17.930 (0.025) 17.441 (0.021) 17.143 (0.014) 16.792 (0.010 ) 

4 16.305 (0.028) 14.861 (0.006) 13.924 (0.016) 12.874 (0.024) 

5 15.624 (0.012) 14.764 (0.028) 14.299 (0.025) 13.850 (0.003) 
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Table 2. Magnitudes of Local Standards: SN 2003du 



Star 



I((Tl 



1 

2 
3 
4 
5 
6 
7 



16.973 
14.944 
17.898 
17.580 
16.431 
13.873 
17.630 
17.977 



(0.013 
(0.025 
(0.022 
(0.032 
(0.026 
(0.025 
(0.028 
(0.031 



16.291 (0.021) 15.875 (0.018 

14.329 (0.017) 13.965 (0.011 

17.025 (0.016) 16.515 (0.015 

16.903 (0.020) 16.508 (0.016 

15.813 (0.021) 15.431 (0.014 

13.351 (0.024) 13.004 (0.025 

16.295 (0.029) 15.339 (0.027 

17.508 (0.028) 17.206 (0.023 



15.508 (0.023) 
13.614 (0.028) 
16.050 (0.027) 
16.171 (0.018) 
15.119 (0.020) 
12.746 (0.030) 
14.226 (0.019) 
16.849 (0.025) 
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Table 3. Photometric Observations of SN 2002bf 



TIT Date^ 






B ((Tn) 
J-' y-' B ) 








I (aA 


2002-02-22 


-10.07 








17.546(0.452)^= 






2002-02-24 


-8.10 








16.759(0.136)^ 






2002-02-28 


-3.95 








16.303(0.105)" 






2002-03-09 


4.95 


16 


488(0.116) 


16.185(0.119) 


16.057(0.116) 


16 


407(0.123) 


2002-03-12 


7.95 


16 


.750(0.048) 


16.285(0.058) 


16.217(0.029) 


16 


929(0.101) 


2002-03-15 


10.93 


16 


.974(0.044) 


16.553(0.136) 


16.421(0.030) 


17 


110(0.096) 


2002-03-19 


14.92 


17 


.401(0.060) 


16.691(0.040) 


16.657(0.033) 


17 


061(0.145) 


2002-03-25 


20.90 


18 


140(0.102) 


16.998(0.039) 


16.697(0.036) 


16 


957(0.076) 


2002-03-29 


24.99 


18 


217(0.571) 


17.019(0.183) 


16.460(0.178) 


16 


804(0.048) 


2002-04-01 


27.84 


18 


667(0.054) 


17.481(0.051) 


16.876(0.077) 


16 


811(0.072) 


2002-04-04 


30.95 


18 


953(0.034) 


17.760(0.029) 


17.090(0.018) 


17 


101(0.090) 


2002-04-07 


33.89 


19 


106(0.064) 


17.856(0.048) 


17.252(0.042) 


17 


255(0.172) 


2002-04-13 


39.86 


19 


.193(0.072) 


18.042(0.179) 


17.451(0.230) 


17 


383(0.414) 


2002-04-22 


48.86 


19 


.352(0.056) 


18.083(0.051) 


17.752(0.152) 


17 


761(0.115) 


2002-04-30 


56.81 


19 


.334(0.079) 


18.400(0.174) 


18.061(0.120) 







^yyyy-mm-dd. 

^Days since estimated date of maximum B brightness, 2002-03-04.37 UT (HJD 
2,452,337.87), taken at the midpoint of the complete series of exposures. 

"Derived from unfiltered image; see text for details. 
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Table 4. Photometric Observations of SN 2003du 



UT Date^ Day'' B (ctb) V (av) R ((Tr) I (ct/) 



2003- 


05- 


01 


-4.74 


















13 


795 


;o 


085) 


13 


938 


;o 


035) 


2003- 


05- 


-06 


0.22 


13 


508 


;o 


028) 


13 


602 


;o 


052) 


13 


572 


;o 


035) 


13 


920 


;o 


031) 


2003- 


05- 


■10 


4.27 


13 


613 


;o 


036) 


13 


637 


;o 


019) 


13 


611 


;o 


015) 


14 


058 


;o 


034) 


2003- 


05- 


■11 


5.25 


13 


650 


;o 


017) 


13 


665 


;o 


019) 


13 


644 


;o 


016) 


14 


105 


;o 


037) 


2003- 


05- 


■12 


6.21 


13 


698 


;o 


033) 


13 


686 


;o 


021) 


13 


674 


;o 


021) 


14 


164 




025) 


2003- 


05- 


13 


7.21 


13 


752 


;o 


031) 


13 


714 


;o 


031) 


13 


729 


;o 


024) 


14 


220 


;o 


034) 


2003- 


05- 


-14 


8.26 


13 


852 


;o 


035) 


13 


759 


;o 


017) 


13 


795 


;o 


033) 


14 


280 


;o 


059) 


2003- 


05- 


15 


9.20 


13 


886 


;o 


031) 


13 


778 


;o 


043) 


13 


847 


;o 


017) 


14 


336 


;o 


028) 


2003- 


05- 


16 


10.23 


13 


981 


;o 


059) 


13 


842 


;o 


024) 


13 


941 


;o 


023) 


14 


427 


;o 


026) 


2003- 


05- 


-17 


11.26 


14 


105 


;o 


131) 


13 


908 


;o 


044) 


14 


004 


;o 


036) 


14 


524 


;o 


028) 


2003- 


05- 


■18 


12.24 


14 


141 


;o 


047) 


13 


973 


;o 


034) 


14 


083 


;o 


020) 


14 


572 


;o 


032) 


2003- 


05- 


■19 


13.27 


14 


246 


;o 


022) 


14 


004 


;o 


019) 


14 


161 


;o 


015) 


14 


638 


;o 


047) 


2003- 


05- 


■20 


14.20 


14 


327 


;o 


016) 


14 


079 


;o 


018) 


14 


227 


;o 


023) 


14 


681 


;o 


044) 


2003- 


05- 


■22 


16.25 


14 


570 


;o 


031) 


14 


197 


;o 


021) 


14 


327 




028) 


14 


710 




032) 


2003- 


05- 


-24 


18.19 


14 


764 


;o 


023) 


14 


298 


;o 


029) 


14 


363 


;o 


026) 


14 


671 


;o 


043) 


2003- 


05- 


26 


20.26 


15 


018 




034) 


14 


404 


;o 


021) 


14 


408 


;o 


020) 


14 


632 


;o 


031) 


2003- 


05- 


28 


22.20 


15 


207 


;o 


016) 


14 


498 


;o 


018) 


14 


428 


;o 


012) 


14 


578 


;o 


033) 


2003- 


05- 


30 


24.23 


15 


404 




039) 


14 


607 


;o 


033) 


14 


439 


;o 


026) 


14 


529 


;o 


031) 


2003- 


05- 


-3lt 


25.22 


15 


565 


;o 


035) 


14 


677 


;o 


014) 


14 


437 


;o 


021) 


14 


467 


;o 


022) 


2003- 


06- 


■01 


26.19 


15 


591 


;o 


027) 


14 


702 


;o 


028) 


14 


475 


;o 


023) 


14 


507 


;o 


038) 


2003- 


06- 


■oit 


26.20 


15 


630 


;o 


032) 


14 


722 


;o 


027) 


14 


468 


;o 


015) 


14 


403 




022) 


2003- 


06- 


■04 


29.21 


15 


866 


;o 


018) 


14 


868 


;o 


024) 


14 


570 


;o 


021) 


14 


454 


;o 


032) 


2003- 


06- 


07 


32.15 


16 


167 


;o 


090) 


15 


054 


;o 


072) 


14 


686 


;o 


030) 


14 


493 


;o 


065) 


2003- 


06- 


10 


35.17 


16 


270 


;o 


025) 


15 


206 


;o 


014) 


14 


883 


;o 


013) 


14 


636 


;o 


037) 


2003- 


06- 


-13 


38.19 


16 


372 


;o 


055) 


15 


388 


;o 


017) 


15 


083 


;o 


014) 


14 


851 


;o 


038) 


2003- 


06- 


■16 


41.13 


16 


511 


;o 


031) 


15 


526 


;o 


034) 


15 


241 


;o 


031) 


15 


036 


;o 


047) 


2003- 


06- 


-19 


44.10 


16 


606 


;o 


015) 


15 


639 


;o 


042) 


15 


374 


;o 


022) 


15 


208 


;o 


051) 


2003- 


06- 


■22 


47.11 


16 


655 


;o 


023) 


15 


748 


;o 


022) 


15 


486 


;o 


015) 


15 


349 


;o 


039) 


2003- 


06- 


■25 


50.09 


16 


705 


;o 


025) 


15 


841 


;o 


022) 


15 


591 


;o 


025) 


15 


511 


;o 


039) 


2003- 


06- 


■26^ 


51.17 


16 


772 


;o 


022) 


15 


842 


;o 


020) 


15 


629 


;o 


037) 


15 


495 


'0 


027) 


2003- 


06- 


■27t 


52.15 


16 


789 


;o 


023) 


15 


872 


;o 


019) 


15 


664 


;o 


027) 


15 


548 




042) 


2003- 


06-28 


53.09 


16 


777 


;o 


025) 


15 


926 


;o 


025) 


15 


697 


;o 


033) 


15 


648 


;o 


045) 


2003- 


07-02 


57.09 


16 


839 


;o 


034) 


16 


031 


;o 


017) 


15 


824 


;o 


021) 


15 


840 


;o 


035) 
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Table 4 — Continued 



UT Date^ 


Day^ 


B (an) 


V (av) 


R {an) 




2003-07-06 


61.09 


16.918(0.017) 


16.135(0.016) 


15.953(0.014) 


16.028(0.040) 


2003-07-10 


65.09 


16.913(0.018) 


16.245(0.017) 


16.088(0.016) 


16.202(0.032) 


2003-07-14 


69.09 


17.052(0.029) 


16.358(0.014) 


16.192(0.022) 


16.331(0.035) 


2003-07-18 


73.08 


17.090(0.024) 


16.456(0.030) 


16.322(0.013) 


16.516(0.030) 


2003-07-22 


77.08 


17.193(0.035) 


16.583(0.019) 


16.461(0.018) 


16.692(0.034) 


2003-08-27t 


113.06 


17.717(0.014) 


17.367(0.054) 


17.506(0.023) 


17.787(0.038) 


2004-07-15* 


435.94 


22.33(0.03) 


22.10(0.04) 


23.44(0.08) 


21.68(0.07) 



Note. — All photometric observations were made with KAIT, except those marked 
with a dagger (^), which were taken with the 1-m Nickel telescope at Lick Observatory, 
and those marked with an asterisk (*), which were taken with HSTs ACS/HRC. 

^yyyy-mm-dd. 

''Days since estimated date of maximum B brightness, 2003-05-06.12 UT (HJD 
2,452,765.62), taken at the midpoint of the complete series of exposures. 



Table 5. Journal of Spectroscopic and Spcctropolarimetric Observations 









HJD 




R angc^' 


Res.'* 


P.A." 


Par. P.A.^ 






See.' 


Slit 


Exp. 




UT Date 


Object 


Day" 


-2,400,000 


Tel.*" 


(A) 


(A) 


(deg) 


(deg) 


Air.8 


Flux Std.*" 


(arcsec) 


(arcsec) 


(s) 


Observers^ 



1997-12-20. 


27 


SN 


1997dt 


20.77: 


50802. 


77 


KIIp 


4310-6820 


5 


98 


77-78 


1.34- 


-1. 


57 


HD84 


1 


.8 


1. 


.0 




4 X 700 


AF, 


AB 




2002-03-07. 


,41 


SN 


2002bf 


3.03 


S2340. 


,91 


Kip 


3900-8800 


13 


98 


172-178 


1.24- 


-1. 


,25 


BD26 


1, 


.5 


1. 


.6 




3 X 300 -H 200 


AF, 


DL, 


EM 


2002-03-11. 


,44 


SN 


2002bf 


7.06 


52344. 


,94 


L3 


3300-6400 


5 


98 


98 




1 


.26 


F34 


2, 


.0 


2. 


.0 




1078 


RC, 


WL 




2002-03-11. 


,44 


SN 


2002bf 


7.06 


52344. 


,94 


L3 


6200-10400 


11 


98 


98 




1, 


.27 


HD19 


2, 


.0 


2. 


.0 




1161 


RC, 


WL 




2003-05-24, 


,46 


SN 


2003du 


18.34 


52783. 


,96 


Kip 


3220-5710 


9 


150 


122-153 


1.35- 


-1, 


,58 


F34 


1 


.5 


1 


,5 


5 


X 900 -1- 3 X 600 


RC, 


DL 




2003-0,5^24. 


46 


SN 


2003du 


18.34 


52783. 


,96 


Kip 


5700-9420 


9 


150 


122-153 


1.35- 


-1, 


,58 


BD17 


1 


.5 


1 


,5 


5 


X 900 -1- 3 X 600 


RC, 


DL 




2003-05-30. 


,40 


SN 


2003du 


24.28 


52789. 


,90 


L3 


3100-5400 


5 


102 


97-101 


1.29- 


-1. 


,32 


F34 


1 


.5 


2. 


.0 




600 + 500 


RF, 


LD, 


DH 


2003-06-30. 


,40 


SN 


2003du 


24.28 


52789. 


,90 


L3 


6200-10400 


11 


102 


97-101 


1.29- 


-1. 


32 


HD84 


1 


.5 


2. 


.0 




600 + 600 


RF, 


LD, 


DH 


2003-06-07. 


39 


SN 


2003du 


32.27 


52797. 


,89 


L3 


3300-6400 


5 


96 


96 




1 


.34 


F34 


2, 


.0 


3 


.0 




700 


RF, 


KS 




2003-06-07. 


39 


SN 


2003du 


32.27 


52797. 


,89 


L3 


6200-10400 


11 


96 


96 




1, 


.34 


BD17 


2, 


.0 


3. 


.0 




700 


RF, 


KS 




2003-07-06. 


,34 


SN 


2003du 


61.22 


52826. 


,84 


L3 


3286-6462 


5 


88 


87 




1 


.46 


F34 


2, 


.0 


2. 


.0 




600 


RF, 


LD, 


MG 


2003-07-06. 


,34 


SN 


2003du 


61.22 


52826. 


,84 


L3 


5150-10200 


11 


88 


87 




1 


.46 


BD26 


2 


.0 


2. 


,0 




600 


RF, 


LD, 


MG 


2003-07-27. 


,20 


SN 


2003du 


82.08 


52847. 


,70 


L3 


3300-5400 


.5 


120 


117 




1 


.18 


BD28 


1 


.0 


2. 


,0 




900 


AF, 


RF, 


LD 


2003-07-27. 


,20 


SN 


2003du 


82.08 


52847. 


,70 


L3 


5200-10400 


11 


120 


117 




1 


.18 


BD17 


1 


.0 


2. 


.0 




900 


AF, 


RF, 


LD 


2004-08-24. 


,43 


SN 


2004dt 


3.93 


53241. 


,93 


L3p 


4600-9800 


18 


160 


145-156 


1.34- 


-1. 


,48 


BD17 


4 


4. 


,0 




4 X 900 


AF, 


RF, 


MG 



^Days since estimated dates of maximum B brightness, for SN 1997dt (1997-11-29; HJD 2,450,782), SN 2002bf {2002-03-04.37; HJD 2,452,337.87), SN 2003du (2003-05-06.12; HJD 2,452,765.62), and 
SN 2004dt (2004-08-20; HJD 2,453,238). The epoch of maximum is particularly uncertain for SN 1997dt, and thus the phase our our single observation is not well known; this uncertainty is indicated by 
a colon. 

^3(p) = Lick 3-m/Kast Double Spectrograph (Miller & Stone 1993); "p" denotes polarimeter attached); KI(II) = Keck-I (II) 10-m/Low-Resolution Imaging Spectrometer (LRIS; Oke et al. 1995); "p" 
denotes polarimeter attached. 

"^Wavelength range of the calibrated flux spectrum. In some cases, the ends are very noisy and are not shown in the figures. 
*^ Approximate spectral resolution derived from night-sky lines. 
^Position angle of the spectrograph slit. 

^Parallactic angle (Filippenko 1982) at midpoint of observation, or PA range for each set of observations. 
^Airmass at midpoint of observation, or airmass range for each set of observations. 

^The standard stars are as follows: F34 = Feige 34, BD28 = BD-^ 28° 4211— Stone (1977), Massey & Gronwall (1990); HD84 = HD 84937, BD26 = BD-K26°2606, BD17 = BD-K17° 4711— Oke &: Gunn 
(1983). 

^Average value of the full width at half maximum of the spatial profile for each set of observations, rounded to the nearest 0^'5. 

■^AB = Aaron Barth; RC = Ryan Chornock; LD = Louis-Benoit Desroches; RF = Ryan Foley; AF = Alex Filippenko; MG = Mohan Ganeshalingam; DH = Deborah Hutchings; WL = Weidong Li; EM 
= Edward Moran; KS = Karin Sandstrom. 
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Table 6. Data on Blueshifts of Si II A6355 



Supernova 


Epoch 


Velocity 
km 


Supernova 


Epoch 


Velocity 
km 


SN 1984A 


-7 


17, 040 




+2 


10, 940 




-6 


16, 840 




+7 


10,400 




-5 


16, 370 




+10 


10,200 




-3 


15, 750 




+14 


10,180 




+8 


13,320 




+21 


10,190 




+17 


12,050 




+28 


9,700 




+ 19 


12,290 


SN 1997bp 


+9 


14, 100 


SN 1991T 


+7 


9,730 


SN 2002bf 


+3 


15,400 




+12 


9,580 




+7 


15, 100 




+13 


9,580 


SN 2002bo 


-12 


17,310 




+18 


9,580 




-1 


14, 000 


SN 1991bg 


-2 


10, 070 




+16 


10, 000 




-1 


9,800 




+28 


9,600 




+17 


7,880 


SN 2003du 


+18 


9,800 




+25 


6,920 




+24 


9,400 


SN 1994D 


-12 


15,580 




+32 


9,500 




-9 


12,530 


SN 2004dt 


-7 


17,200 




-6 


11,100 




+4 


13,500 




-3 


11,160 









Note. — Velocities of absorption minima of the Si II A6355 hne, 
derived using the relativistic Dopplcr formula. All velocities were mea- 
sured directly from spectra in our database, except for SN 1984A, for 
which the Si II minima wavelengths reported by Barbon et al. (1989) 
were used, and the day —7 data for SN 2004dt, for which the wavelength 
of minimum reported by Wang et al. (2005a) was used. All spectra were 
first deredshifted according to the values reported in NED. The adopted 
recession velocities and dates of B maximum are as follows: SN 1984A 
— -yrec = 261 km s~^, B^^ on 1984 Jan. 17 (Barbon et al. 1989); 
SN 1991T — v,ec = 1736 km s'^, B^^ on 1991 Apr. 28 (Phillips et al. 
1992); SN 1991bg — v,^ = 1060 km s-^ B^^ on 1991 Dec. 15 (Lei- 
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bundgut et al. 1993); SN 1994D — v^ec = 850 km s'K B,,,^^ on 1994 
Mar. 20 (Richmond et al. 1995); SN 1997bp — v,^^ = 2492 km s-\ 
B^^ on 1997 Apr. 7 (Riess et al. 1998b); SN 2002bf — i;rec = 7254 km 
Smax on 2002 Mar. 4 (this work); SN 2002bo — v,^ = 1271 km s-^ 
B^^ on 2002 Mar. 23 (Benetti et al. 2004); SN 2003du — v^ec = 1194 
km s"\ Smax on 2003 May 6 (this work); SN 2004dt — -yrec = 5915 km 
s~^, -Bmax on 2004 Aug. 20 (W. Li, private communication). 
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Table 7. MLCS2k2 Parameters of SN 2002bf and SN 2003du 











My (a)'' 


A^(a)d 




B f 

-'-'max 


V(at Bmax)f 


SN 


(mag) 


(Mpc) 


A(a)b 


(mag) 


(mag) 


Ma)" 


(mag) 


(mag) 


SN 2002bf 


35.47(0.13) 


124(8) 


-0.11(0.07) 


-19.58(0.05) 


0.26(0.13) 


52337.87(0.50) 


16.22 


16.18 


SN 2003du 


33.18(0.05) 


43(1) 


-0.25(0.03) 


-19.67(0.02) 


0.04(0.03) 


52765.62(0.50) 


13.56 


13.58 



^Statistical uncertainty only. 

''MLCS2k2 light-curve shape parameter, with lower values of A corresponding to brighter, more slowly declining light curves, 
relative to the "fiducial" A = template, which has Ami5{B) = 1.07 mag. 

"^y-band absolute magnitude at the time of peak B-band brightness, on the Hq = 65 km s ~^ Mpc~^ distance scale used 
here. The fiducial A = template has the following peaJc magnitudes at B maximum light: Mg = —19.56 mag. My = —19.50 
mag, = -19.50 mag, and Mf = -19.20 mag. 

'^Host-galaxy extinction, assuming an Ry = 3.1 extinction law. 

''Heliocentric Julian date minus 2,400,000 for the time of B-band peak. The time of y-band maximum for both SNe is about 

1 day after the time of B-band peak. 

^Observer frame, no fe-correction, no extinction correction. V-band peak at t = t(ymax) is 0.01 mag brighter than y(at to) 
for both objects. 



